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SUMMARY 

It  has  been  established  from  earlier  studies  that  aircraft  engine 
pollutant  emissions  are  becoming  a  more  significant  factor  in  regional 
pollutant  le.'els  and  are  already  of  concern  in  the  communities  adjacent 
to  the  major  air  terminals.  The  pollutant  classes  to  be  contended  with 
have  been  shown  to  be  carbon  monoxide,  the  hydrocarbons,  nitrogen  ox¬ 
ides  (NOx),  and  particulates.  Of  these,  the  only  pollutant  class  for 
whi ch  .control  criteria  do  not  exist  are  the  nitrogen  oxides,  and  the 
object  of  this  study  was  to  produce  such  criteria  as  might  apply  to  the 
range  of  aircraft  gas  turbine  engines  that  are  in  use,  or  i~  the  plan¬ 
ning  phase  today. 

The  method  of  apprcacn  adopted  to  develop  the  des:gn  criteria  in¬ 
volved  four  separate  studies  with  the  following  objectives: 

a.  The  development  of  a  computerized  mathematical  model  of  NOx 
emission  from  aircraft  engine  combustors. 

b.  A  parametric  analysis;  using  the  model,  to  determine  the 
sensitivity  of  NOx  emission  to  variations  oc  model  parameters 
and  engine  design  variables. 

c.  Evaluation  of  critical  model  parameters  by  means  of  experi¬ 
mental  measurements. 

d.  The  incorporation  of  the  model  into  existing  combustor  design 
methods  to  provide  guidelines  for  minimizing  NOx  emission 
while  maintaining  other  performance  characteristics  at  satis¬ 
factory  levels 

The  mathematical  model  presented  is  based  upon  a  relatively  simple 
approach  to  the  problem  and  its  validity  rests  upon  two  main  assumptions 
firstly,  that  the  combustion  process  is  mixing  controlled  and  therefo'e, 
that  equilibrium  conditions  occur  for  the  main  species  or  interest  in 
the  nitric  oxide  formation  process;  and  secondlv,  that  the  degre.  of  i 
perfect  mixing  in  the  hot  zones  of  the  combustor  can  be  represented 
statistically  by  a  Gaussian  distribution  in  fuel-air  ratio.  Based  on 
these  assumptions,  the  model  predicts  distribution  character i si i cs  and 
the  mass  average  values  of  nitric  oxid.  concentrations,  temperature. 
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fuel  concent rat  ions,  and  residence  times  at  all  regions  throughout  the 
combustor. 

The  computer  model  was  first  tested  with  a  parametric  study  and 
it  involved  both  a  model  sensitivity  analysis,  which  established  the 
relative  importance  of  certain  parameters  such  as  the  mixing  parameters 
which  are  internal  to  the  computing  procedure,  and  then  a  design 
sensitivity  analysis.  This  latter  analysis  chiefly  established  the 
importance  of  the  combustor  operating  conditions  and  showed  that  the 
engine  compressor  pressure  ratio  (or  combustor  i nleti  temperature)  has 
the  greatest  influence  upon  the  nitric  oxide  emission  levels  for  a 
given  a:rcraft  operating  condition.  Of  the  variables  under  the  control 
of  the  designer,  the  mass  mean  equivalence  ratio  and  the  degree  of  homo¬ 
geneity  of  combustion  within  the  hot  primary  zone  of  the  combustor  were 
established  tc  be  the  most  significant  to  control,  with  minimum  emis¬ 
sion  levels  corresponding  to  the  lowest  equivalence  ratio  and  best 
mixed  condition. 

The  predictions  of  the  computer  model  were  then  tested  against 
new  experimental  data  collected  during  this  study.  The  results  pre¬ 
sented  are  for  two  current  aircraft  engine  combust,  rs,  designed  by 
different  design  teams,  that  were  operated  at  conditions  which  simulated 
those  experienced  in  flight  operation.  The  conditions  applied  to  each 
combustor  were  made  to  overlap  so  that  their  nitric  oxide  emission  charac 
teristics  could  be  compared.  Significant  differences  can  be  observed 
between  these  combustors,  and  the  results  are  found  tr  correlate  very 
well  with  the  model  predictions  for  certain  prescrioed  values  of  the  homo 
geneity  of  combustion  in  the  primary  zone.  The  homogeneity  parameter  is 
found  not  to  dominate  the  predictions  and  the  predicted  results  are  found 
to  be  influenced  more  by  the  other  combustor  operating  conditions. 

Based  on  the  results  of  these  studies  and  predictions  made  by  the 
mathema* i ca!  model,  a  design  technique  is  established  which  may  be  used 
to  minimize  future  a:rcraft  gas  turbine  engine  emissions  of  nitric  oxide. 
Gap  areas  in  know'edge  about  combustor  operatirg  phenomena  currently 
limit  the  usefulness  of  the  technique,  but  it  has  been  applied  to  one 
of  the  combustors  tested  to  estimate  the  potential  reduction  in  nitric 
oxide  emissions  that  might  be  attained  if  the  technique  were  employed. 
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Alternative  methods  c'.  nitric  oxide  control  are  considered  and  the 
potential  effectiveness  of  each  control  method  is  evaluated  in  certain 
cases.  The  alternative  control  methods  considered  are  water  injection, 
variable  geometry,  and  staged  combustion. 

On  the  basis  of  the  results  of  this  study,  the  following  conclu¬ 
sions  have  been  formulated: 

a.  NGx  emission  rate  can  be  related  to  fundamental  combustor 
design  parameters  which  also  influence  other  major  combustion 
performance  factors-  This  interrelationship  between  NOx 
emissions  and  other  combustion  performance  factors  provides 

a  basis  for  quantitative  NOx  emission  control  criteria. 

b.  NOx  emission  rates  from  conventional  combustors  can  be  reduced 
by  modification  of  fuel  distribution  parameters.  However, 

NOx  reduction  may  be  accompanied  by  significant  reductions 

in  other  performance  factors  such  as  combustion  efficiency 
and  reignition  altitude.  These  performance  losses  are  more 
severe  in  a  mod if :cd-ccmbustor  design  procedure  where  com¬ 
bustor  size  is  fixed  than  in  a  new-combus tor  design  procedure. 

c.  Water  injection  =  s  effective  in  reducing  NOx  emissions  from 
conventional  combustors.  Emission  rate  reductions  of  approxi¬ 
mately  50  per  cent  at  high  power  levels  can  be  obtained  with 
water  injection  rates  equal  to  fuel  injection  rates.  Increased 
control  may  be  obtained  with  increased  water  injection  rates, 
or  by  localizing  injection  of  water  into  the  combustion  zone. 

d.  Advanced  combustor  design  concepts,  such  as  variable  geometry 
or  staged  combustion,  show  .oromise  of  effective  NOx  emission 
control  for  future  engines.  The  degree  of  NOx  emission  control 
attainable  using  either  of  these  concepts  is  estimated  to  be 
of  the  same  order  as  that  attained  by  water  injection. 

A’l  the  results  of  this  study  are  presented  in  Volumes  1,  2,  and 
3  of  this  report.  Volume  1  describes  the  complete  program,  but  precise 
derails  o*  the  mathematical  model  and  the  computer  programs  developed 
from  the  model  are  presented  in  the  other  volumes.  Volume  2  describes 
the  theory  and  computer  program  for  determining  nitric  oxide  formation 
rates,  and  Volume  3  the  corresponding  information  for  the  flow  model  used 
to  determine  the  nitric  oxide  emission  level. 
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The  rates  of  pollutant  emission  by  aircraft  are  increasing  with 
time  as  a  result  of  increases  in  aircraft  activity  and  certain  changes 
in  equipment,  including  size,  engine  type,  and  performance  characteris¬ 
tics.  Thus,  aircraft  emissions  are  becoming  a  more  important  factor 
in  regional  pollution  levels.  This  trend  will  continue  until  criteria 
are  developed  for  reducing  emissions,  and  standards  are  set  to  limit 
emission  rates. 

Design  methods  exist  today  for  minimizing  emissions  of.  carbon 
monoxide  and  hydrocarbons  (Ref  4).  The  successful  application  of  these 
methods  at  medium  and  high  power  operation  conditions  has  resulted  in 
high  efficiency  and  low  emission  rates  during  aircraft  flight  operations. 
The  application  of  these  design  methods  to  low  power  conditions  should 
reduce  emissions  of  CO,  hydrocarbons,  and  possibly  odor  by  aircraft 
engines  during  ground  operations.  Methods  of  reducing  these  "low 
power  emissions"  are  being  developed  currently  by  the  Pratt  &  Whitney 
Aircraft  Division  of  the  United  Aircraft  Corporation  for  the  Air 
Force  Aero  Propulsion  Laboratory. 

Design  methods  of  an  empirical  nature  exist  for  reducing  the 
emission  of  particulates  by  aircraft  turbine  engines.  These  methods 
have  been  utilized  in  the  past  in  the  development  of  combustion 
chambers  in  order  to  achieve  maximum  chamber  durability  and  to  reduce 
carbon  deposition  throughout  the  engine-  More  stringent  applications 
of  these  same  methods  are  serving  to  reduce  visible  smoke  emission  from 
contemporary  engines. 

The  only  pollutants  for  which  control  criteria  do  not  exist  are 
the  nitrogen  oxides.  Because  of  this  lack  of  NOx  control  criteria, 
it  is  not  possible  to  predict  accurately  the  effects  o*  engine  design 
changes  on  rates  of  NOx  emission.  However,  it  is  clear  that  current 
trends  toward  increased  engine  power  levels  and  operating  temperatures 
will  result  in  higher  NOx  emission  rates  unless  control  criteria  are 
■developed. 

2.2  PROGRAM  OBJECTIVE 

The  objective  of  this  study  was  to  develop  criteria  which  may  be 
used  during  the  design  of  aircraft  gas  turbine  combustion  chambers  to 
minimize  the  emissions  of  nitrogen  oxides. 


To  be  useful,  design  criteria  for  NOx  control  must  provide 
quantitative  relationships  between  design  variables  and  emission  rates. 
Such  criteria  can  be  developed  by  empirical  analysis  of  large  quantities 
of  engine  design  and  emission  data,  or  by  theoretical  analysis  of  com¬ 
bustor  performance  supported  by  limited  experimental  data.  The  method 
adopted  for  this  program  is  based  on  the  latter  approach,  since  it  is 
considered  to  be  a  more  effective  and  economical  approach  to  this 
particular  problem. 

Having  selected  the  theoretical  analysis  approach,  two  alterna¬ 
tives  present  themselves-  A  highly  accurate  and  comprehensive  model 
of  NOx  emission  might  be  selected  as  the  program  goal,  in  which  case 
complete  achievement  of  the  goal  in  the  initial  program  could  not  be 
assured.  Alternatively,  one  might  set  out  to  deveJop  a  simpler  but 
still  valid  NOx  emission  model  and  then  to  evolve  design  criteria 
based  upon  this  model.  This  latter  alternative  was  selected  as  a  goal 
which  could  be  achieved  with  a  reasonable  amount  of  effort. 

Thus,  the  problem  to  which  this  study  was  addressed  is  the  devel¬ 
opment  of  criteria  for  relating  NOx  emissions  to  engine  design  and  per¬ 
formance  characteristics,  and  the  method  of  approach  adopted  involved 
four  separate  studies  with  the  following  objectives: 

a.  The  development  of  a  computerized  mathematical  model  of  NOx 
emission  from  aircraft  engine  combustors. 

b.  A  parametric  analysis,  using  the  model,  to  determine  the 
sensitivity  of  NOx  emission  to  variations  of  model  parameters 
and  engine  design  variables. 

c.  Evaluation  of  critical  model  parameters  by  means  of  experi¬ 
mental  measurements. 

d.  The  incorporation  of  the  model  into  existing  combustor 
design  methods  to  provide  guidelines  for  minimizing  NOx 
emission  while  maintaining  other  performance  and  emission 
characteristics. 

Volume  1  of  this  report  describes  the  methods  used  to  satisfy 
these  objectives  and  also  the  results  obtained  from  all  four  tasks 
described  above.  Details  of  the  mathematical  model  and  the  computer 
program  developed  from  the  model  are  fully  described  in  Volumes  2 
and  3* 
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3.  THE  EMISSIONS  MODEL 

Calculations  based  upon  NOx  kinetics  clearly  show  that  for  mea¬ 
surable  concentrations  of  nitric  oxide  to  be  formed  during  the  short 
time  it  takes  the  gases  to  pass  through  an  aircraft  gas  turbine  engine, 
the  temperature  must  exceed  2000  deg  K.  Such  high  temperatures  only 
exrst  within  the  combustor,  and  only  then  for  a  limited  time,  so 
clearly,  nitric  oxide  concentrations  in  the  exhausts  of  aircraft 
engines  are  solely  dependent  upon  the  flow  behavior  and  the  chemical 
processes  occurring  within  the  combustion  chamber  and  it  is  these 
features  that  need  to  be  adequately  represented  in  the  model-  The 
important  chemical  reactions  and  their  reaction  rates  must  be  known 
and  a  rate  equation  for  nitric  oxide  formation  derived  from  these 
data.  A  flow  model  also  has  to  be  developed  to  predict  the  thermo¬ 
dynamic  states  of  the  gases  inside  the  combustor  liner.  The  gas  tem¬ 
perature,  density,  and  concentrations  of  species  important  in  the 
nitric  oxide  formation  process  must  then  be  determined  and  the  rate 
equation  integrated  through  this  sequence  of  thermodynamic  states. 

Conveniently,  the  two  parts  can  be  separated  for  two  reasons- 
Firstly,  nitric  oxide  concentrations  are  sufficiently  low  so  that 
changes  in  concentration  have  a  negligible  effect  upon  the  gas  tem¬ 
perature.  Secondly,  the  nitric  oxide  formation  process  proceeds  more 
slowly  than  fuel  combustion  processes  so  that  most  nitric  oxide  ir 
formed  after  concentrations  of  other  combustion  products  have  reached 
equilibrium  levels-  The  nitric  oxide  formation  process  and  combustor 
flow  processes  are  considered  separately  below. 

3-1  THE  NITRIC  OXIDE  FORMATION  PROCESS 

Engine  exhaust  emissions  are  often  expressed  as  nitrogen  oxides 
(NCx) ,  or  as  nitrogen  dioxide  (NO^)-  However,  at  the  high  power 
conditions,  the  pollutant  present  in  the  exhaust  gas  as  it  leaves  the 
engine  is  almost  completely  nitric  oxide  (NO,  with  only  a  few  per  cent 
of  NO^  (Ref  5)*  At  low  engine  power  levels,  the  percentage  of  N0^ 
present  can  increase  significantly  but  at  these  conditions  the  total 
NOx  concentration  is  very  low,  hence,  the  NO^  contribution  is  not 


significant  when  compared  to  levels  of  total  HOx  emissions  which  are 
of  interest  today-  The  conversion  to  NOx  or  NO^  is  made  because  many 
of  the  early  techniques  used  to  measure  nitric  oxide  required  first 
that  it  be  oxidized  to  nitrogen  dioxide.  Thus,  inside  the  engine, 
the  important  oxide  of  nitrogen  to  consider  is  nitric  oxide,  NO. 

There  are  many  reactions  schemes  proposed  in  the  .iterature  for 
nitric  oxide  formation  (Refs  5-8)  and  they  are  reviewed  in  Volume  2, 
Section  2  of  this  report  and  in  Reference  9-  The  reaction  scheme 
employed  in  these  studies  and  incorporated  into  the  flow  model  is  as 


shown  below: 

N  +  NO  ^  N2  -r  0  (1) 

N  +  02  — NO  +  0  (2) 

N  +  OH  v  NO  +  H  (3) 

ri  +  N20  L - T  N2  +  OH  (*0 

0  +  n2o^=?  N2  +  o2  (5) 

0  +  N20  £=y  NO  +  NO  (6) 


Each  of  the  six  reactions  proceeds  at  »  different  reaction  rate 
which  can  be  defined  by  a  general  rate  term  of  the  form 

k.  =  A.Tn'  exp (-E./RT)  (7) 

where  k.  is  the  reaction  rate  constant  for  the  eauatior.  i 

i 

T  is  temperature 
R  is  the  universal  gas  constant 

and  A,  n,  and  E  are  empirical  constants  which  need  to  be  determined 
by  experiment.  Many  values  are  given  for  these  constants  in  the 
literature  and  these  too  are  reviewed  in  Volume  2,  Section  2,  and 
four  different  sets  of  data  are  isolated. 

In  general,  one  can  show  that  not  all  of  the  above  reactions 
contribute  significantly  to  the  formation  of  nitric  oxide  over  the 
conditions  encountered  in  the  aircraft  gas  turbine  engine.  The 
contributions  by  the  first  two  reactions  (Equations  1  and  2),  often 
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referred  to  as  the  ZeJdovich  chain  mechanism,  are  by  far  the 
most  significant.  The  third  reaction  can  become  important  under 
particular  operating  conditions,  such  as  when  parts  of  the  fuel 
charge  burn  in  the  presence  of  an  insufficient  quantity  of 
air,  but  simple  calculations  show  that  reactions  4,  5,  snd  o 
do  not  contribute  greatly  to  the  nitric  oxide  formation  rate. 
However;  in  the  interest  of  completeness,  all  reactions  have 
been  included  in  the  computerized  model  and  their  individual 
contributions  studied  as  part  of  the  second  phase  of  this  program 
(see  Section  k. 1). 

In  order  to  apply  this  reaction  scheme  to  gas  turbine 
combustors,  it  is  necessary  to  be  able  to  predict  the  concen¬ 
tration  levels  of  all  species  active  in  the  set  of  six  reactions 
specified  above  at  all  thermodynamic  states  typical  of  the  com¬ 
bustor.  At  the  pressure  and  temperature  levels  that  exist  in  such 
combustors  during  high  thrust,  low  altitude  operation,  the  hydro¬ 
carbon  oxidation  reactions  go  rapidly  to  completion  compared  with 
the  nitric  oxide  formation  processes.  In  order  to  predict  these 
specie  concentrations  therefore,  the  following  assumptions  are 
made: 

a.  Combustion  is  mixing  controlled  and  not  reaction  rate- 
limited,  and  therefore,  that  in  the  combustor,  N2,  0 2, 

0,  OH,  H,  and  H2  concentrations  are  the  equilibrium  con¬ 
centrations  corresponding  to  the  combustor  inlet  tempera¬ 
ture  and  pressure,  and  to  the  equivalence  ratio.  This 
assumption  limits  the  usefulness  of  the  model  to  high 
pressure,  high  temperature  operating  conditions,  but  it 
is  precisely  these  conditions  at  which  most  nitric  oxide 
is  produced  and  are  therefore  of  most  interest  to  this 
study. 

b.  N  and  N20  concentrations  are  given  by  the  steady-state 
assumption;  that  these  concentrations  are  not  in  equilib¬ 
rium  but  are  in  steady-state  with  NO  can  be  shown  by 
deriving  rate  equations  for  N  and  N20  which  have  relax¬ 
ation  times  short  compared  with  the  NO  rate  equation. 


If  these  assumptions  are  applied  to  Equations  ■  through  7i  then 
the  change  in  total  nitric  oxide  per  unit  time  per  unit  volume  that 
results  from  chemical  reaction  over  a  fixed  element  can  be  expressed 
as  a  function  solely  of  the  local  nitric  oxide  concentration  and  the 
local  pressure,  temperature,  and  chemical  specie  composi t ions  (see 
Equation  8,  Volume  2,  Section  2).  If  a  flow  model  is  now  developed 
which  can  predict  values  of  these  latter  variables  throughout  the  gas 
turbine  engine,  then  the  generation  of  nitric  oxide  can  be  calculated 
by  integration  of  the  chemical  rate  function. 

THE  FLOW  MODEL 

The  flow  model  under  consideration  has  to  he  applied  to  both 
canannular  and  annular  combustors  and  to  take  into  account  the  two 
distinct  functions  performed  in  their  operation.  These  are  the  rapid 
mixing  between  fuel  and  air,  flame  stabilization,  and  combustion  pro¬ 
cesses  occurring  in  the  primary  zone,  and  the  combustion  gas  cooling 
by  mixing  with  dilution  air  before  entry  into  the  turbine  in  the 
secondary  zones.  Figure  la  is  a  schematic  of  the  maan  flow  pattern 
illustrating  these  two  functions.  The  flow  characteristics  in  each 
zone  are  significantly  different. 

In  the  primary  zone,  a  mean  flow  pattern  with  recirculation  approxi¬ 
mately  as  shown  in  Figure  la  can  be  defined.  However,  substantial  tur¬ 
bulent  fluctuations  about  this  mean  occur,  resulting  in  wide  variations 
in  residence  time  and  composition  within  the  zone-  The  fuel,  air,  and 
combustion  products  are  not  uniformly  mixed,  and  pBrts  of  the  flow  are 
significantly  leaner  and  parts  are  richer  than  the  mean  fuel -air  ratio. 
One  result  of  this  incomplete  mixing  is  that  some  unreacted  fue!  leaves 
the  primary  zone.  Even  at  full  power,  low  altitude,  operation  when  *he 
kinetics  of  the  hydrocarbon  oxidation  process  are  not  important,  less 
than  90  per  cent  of  the  fuel  fed  to  the  primary  zone  is  burned  there. 


The  term  "burn"  is  used  throughout  this  report  to  represent  the 
processes  of  ignition  and  reaction  of  fuel -air  mixtures  to  form  com¬ 
bustion  products  of  equilibrium  composition  (except  for  NO).  Burning 
is  assumed  to  occur  instantaneously  upon  mixing  of  fuel  and  air 
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In  the  rest  of  the  combustor,  the  remaining  fuel  is  burned,  and 
diluting  air  is  mixed  with  the  primary  zone  products-  While  turbulent 
fluctuations  are  significant,  the  flow  is  primarily  un  directional. 
Typically,  dilution  air  is  added  both  through  large  holes  in  the  liner 
so  that  air  jets  penetrate  to  the  core  of  the  combustion  products  flow, 
through  smaller  holes  so  that  mixing  occurs  at  intermed:ate  radii,  and 
along  the  walls  for  film  cooling-  Thus,  at  combustor  exit,  the  gases 
are  sufficiently  well  mixed  to  obtain  the  desired  temperature  profile 
across  the  flow  before  entry  to  the  turbine.  Separate  flow  models 
are  required  for  these  different  parts  of  the  combustor.  Figure  1 
illustrates  the  division  of  the  liner  into  three  zones  and  the  model 
used  in  each  zone.  These  models  are  described  qualitatively  below 
and  full  mathematical  detaiis  are  presented  in  Volume  3»  Section  2,  of 
this  report. 

3-2-1  Primary  Zone  Model 

-The  hottest  zone  in  a  gas  turbine  combustor  is  normally  the 
primary  zone,  it  is  the  zone  where  fuel  and  air  are  injected  at  an 
approximately  stoichiometric  mixture  strength  and  where  most  of  the 
combustion  takes  place.  The  fuel  is  sprayed  from  an  atomizing  fuei 
nozzle  and  strategically  located  jets  cause  the  incoming  air  to 
generate  a  recirculating  flow  of  the  type  shown  schematically  in 
Figure  la.  It  is  a  region  of  intense  turbulent  combustion  and  the 
process  is  so  complex  that  it  is  not  possible  to  produce  a  ccxnplete 
analytical  (i-e-,  mathemat ical )  description  of  the  time  sequence  of 
events. 

One  analytical  approach  used  to  represent  reacting  systems  such 
as  that  above  is  to  adopt  the  concept  of  a  perfectly  stirred  reactor. 
That  is  to  assume  that  the  reactants  mix  instantaneously  to  the 
perfectly  mixed  state  before  reaction  proceeds.  It  is  clear  that 
such  a  reactor  is  an  ideal  concept  that  can  only  be  approached  in 
practical  systems  since  mixing  develops  at  a  finite  rate.  The  appli¬ 
cation  of  this  concept  to  the  problem  under  study  can  be  anticipated 
to  be  totally  inadequate  since: 
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a.  The  mixing  orocess,  involving  reactants  in  different  phases, 
takes  a  long  time  relative  to  the  combustion  process  and 
causes  combustion  to  take  place  between  imperfectly  mixed 
reactants. 

b.  Primary  zone  equivalence  ratios  are  typically  about  unity. 
Distributions  in  the  fuel-air  ratio  about  this  mean  value 
substantially  change  the  concentrations  of  the  major  chemical 
species  and  temperatures,  and  can  be  shown  (from  Equations  1 
through  7)  to  drastically  influence  nitric  oxide  formation 
rates. 

Hence,  the  effect  of  imperfect  mixing  must  be  taken  properly  into 
account. 

In  addition  to  *he  calculation  of  nitric  oxide  formation  rate,  it 
is  also  necessary  to  determine  the  time  the  reactants  spend  in  the  pri¬ 
mary  zone,  that  is  the  residence  time.  For  the  perfectly  stirred 
reactor,  this  presents  no  problem  as  one  can  show  mathematically  that 
there  must  be  a  distribution  of  residence  times  in  an  exponential  form 
related  to  the  volumetric  flow  rate  and  the  volume  of  the  reactor  (see 
Volume  3,  Section  2).  This  same  distribution  of  residence  times  will 
be  assumed  to  be  applicable  as  it  has  been  shown  that  most  practical 
systems  do  in  fact  exhibit  this  property. 

The  problem  in  developing  an  adequate  model  for  the  combustion 
in  the  primary  zone,  therefore  becomes  one  of  accounting  for  the 
nonuniformities  in  concentration  that  exist  at  the  onset  of  combustion. 

It  is  assumed  here  that  the  flow  into  the  primary  zone  rapidly  mixes 
to  a  scale  small  compared  wi th  combustor  dimensions,  but  not  complete 
on  a  molecular  scale,  such  that  there  exist  discrete  eddies  of  different 
mixtures  which  burn  and  retain  their  identity  throughout  their  passage 
through  the  primary  zone.  As  the  scale  of  mixing  is  assumed  smal1,  it 
is  consistent  to  assume  that  the  nonuniformi *■  ies  in  concentration  can 
be  represented  by  a  Gaussian  distribution  function  as  shown  in  Figure  lc- 

The  "mixedness"  of  the  primary  zone  can  then  be  characterized  by 
the  single  dimensionless  parameter  Sq  =  /  pj,  where  CT  in  statis¬ 

tical  terms  is  called  the  standard  deviation  and  represents  the  degree 
of  distribution  of  F  about  F  so  that  the  quality  of  mixing  decreases 


as  the  parameter  S  increases,  (it  should  be  noted  that  S  =0 
°  o 

corresponds  to  the  perfectly  mixed  case,  i.e. ,  to  the  perfectly 

stirred  reactor.)  At  the  present,  the  value  of  S  can  on’y  be  esti- 

o 

mated,  bet  clearly  it  ?s  related  to  the  energy  supplied  to  the 
primary  zone  by  the  air  and  fuel  jets.  Its  value  together  with  those 
of  other  design  variables,  is  first  treated  parametr ical ly  in  the 
initial  computer  analysis  presented  in  Section  h,  and  one  of  the  ob¬ 
jectives  of  the  experimental  program  described  in  Section  5  is  to 
isolate  its  value  for  the  two  combustor  tests. 

As  previously  discussed,  not  sli  the  fuel  supplied  to  the  primary 
zone  actually  burns  in  the  zone.  Account  is  taken  of  this  fact  simply 

by  modifying  the  mass-mean  mixture  ratio  F  .  calculated  from  the  fuel 

P* 

and  airflow  rates  into  the  zene.  The  fraction  of  fuel  burned,  repre¬ 
sented  by  the  symbol  j3  ,  is  calculated  from  the  correlation  shown  in 
Figure  9,  taken  from  Reference  i5-  This  correlation  is  based  upon  mea¬ 
sured  primary  zone  combustion  efficiency  (T?  ).  Since  combustion 

'  cp 

efficiency  is  reduced  both  by  nonburning  and  partial-turning  of  fuel, 
the  use  of  this  correlation  to  calculate  the  fraction  of  fuel  burned 
introduces  an  error  into  the  analysis.  The  effects  of  partial-burning 
of  fuel  are  introduced  through  the  mixture  ratio  distribution  factors 
Fp  and  Sq  since  partial-burning  results  primarily  from  the  use  of 
rich  mixtures  or  poor  mixing.  Therefore, 0  should  be  evaluated  on 
the  basis  of  nonburnsd  fuel  only.  However,  data  on  nonburned  fuel  at 
the  primary  zone  exit:  as  distinguished  from  primary  zone  efficiency, 
are  not  available.  Hence,  the  efficiency  correlation  'S  used  to  cal¬ 
culated  and  the  error  thus  introduced  is  noted;  that  is,  qualitatively, 
&  should  be  somewhat  higher  than  '»)  .  Tne  unburned  fuel  is  assumed 

*  CO 

to  burr,  in  the  intermediate  zone  in  a  manner  as  described  be’ow. 

Combustion  dimensions,  airflow  distribut Son  characteristics,  oper¬ 
ating  conditions,  and  a  prescribed  value  of  the  mixedness  parameter  are 
all  that  are  needed  to  calculate  the  mass  average  nitric  oxide  concen¬ 
tration  at  the  exit  from  the  primary  zone. 
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3-2.2  The  Intermediate  Zone 

The  intermediate  zone  represents  a  transition  between  the  primary 
and  dilution  zones  as  shown  in  Figure  la-  Mixing  occurs  in  this  zone 
between  the  heterogeneous  products  from  the  primary  zone  and  the 
entering  secondary  combust i on  air  which,  in  practice,  is  usually  added 
sufficiently  slowly  to  allow  recombi nation  of  the  dissociated  combustion 
products  and  burn-up  of  the  remaining  unburned  fuel  to  occur.  The  basic 
assumption  made  in  modeling  the  flow  is  that  In  the  prescribed  length, 

L,  complete  mixing  occurs  between  ail  elements  leaving  the  primary  zone, 
secondary  air  added  to  the  zone,  and  all  fuel  not  burned  in  the  primary 
zone.  Entry  conditions  to  this  zone  correspond  to  time>-average  primary 
zone  exit  conditions;  at  zone  exit  the  flow  is  homogeneous  and  the  flow 
velocity  is  assumed  to  be  one-dimens  tonal .  it  should  be  noted  that 
the  value  of  L  can  be  greater  than  the  length  of  the  combustor,  in 
which  case  complete  mixing  will  not  occur  within  the  system. 

The  mixing  conditions  in  this  zone,  can  be  represented  schemati¬ 
cally  as  shown  in  Figure  lb.  The  mass-average  mixture  ratio  usually 
decreases  from  that  at  the  primary  zone  exit  due  to  the  addition  of 
moo  air  through  holes  in  the  liner  wail  and  the  degree  of  unmixedness 
steadily  decreases.  In  order  to  follow  the  change  in  nitric  oxide  con¬ 
centration  along  the  zone,  the  following  assumptions  are  mace  concerning 
the  mixing  behavior  (see  Volume  3,  Section  2): 

a.  That  the  rate  of  mixing  to  the  homogeneous  state  is  propor¬ 
tional  to  the  square  root  of  the  distance  traveled  from  the 
primary  zone  exit. 

b.  That  the  unburned  fuel  leaving  the  primary  zone  follows  the 
same  mixing  law  as  for  (a)  and  is  instantaneously  burned  at 
the  time  of  mixing. 

c.  That  air  entering  the  intermediate  zone  through  the  liner 
wall  is  distributed  among  ail  fiow  elements  in  proportion  to 
the  masses  of  the  elements,  and  that  s.ixing  of  this  air  occurs 
at  a  finite  rate  according  to  the  same  mixing  law  as  for  (a). 

d.  That  intermixing  occurs  between  elements  of  different  nitric 
oxide  concentration,  and  the  rate  of  mass  flow  from  an  element 
due  to  intermixing  is  proportional  to  the  mass  of  the  element- 
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Justification  for  the  first  two  assumptions  rests  with  the  fact 
that  the  mixing  rate  is  known  to  be  proportional  to  /T,  i.e.  ,  is 
approximately  proportional  to  '/x  and  it  is  this  type  of  mixing  which 
is  finally  of  importance  in  the  reacting  process.  The  latter  two 
assumptions  stem  from  the  concept  of  combustor  aerodynamics  employed 
in  this  study  in  which  the  flow  system  is  envisaged  as  a  multitude 
of  eddies  each  of  which  retains  its  identity  until  mixing  is  complete. 
All  eddies  are  assumed  to  be  similar  in  size  (or  mass)  and  to  exchange 
mass  with  other  eddies  through  mixing  processes  at  similar  rates.  A 
flow  element  includes  all  eddies  with  fuel-air  ratios  within  a  specified 
range.  Hence,  it  follows  that  the  rate  of  exchange  of  mass  among  ele¬ 
ments  is  proportional  to  the  masses  of  the  elements.  The  sensitivity 
of  the  results  to  this  assumed  mixing  behavior  is  explored  In  Section  4 
of  this  report. 

3*2-3  The  Dilution  Zone 

in  *he  dilution  zone,  the  flow  is  assumed  to  be  one-dimensional 
with  tt  gases  uniformly  mixed  across  each  cross  section.  Only  the 
mean  mixture  ratio  need  be  considered,  and  this  only  changes  as  the 
remaining  compressor  air  is  mixed  in.  The  nitric  oxide  mass  fraction 
changes  through  chemical  reaction,  and  through  dilution  as  additional 
air  ! s  mixed  with  the  bulk  flow  within  the  combustor  liner  and  both 
factors  can  be  taken  properly  into  account  by  the  method  presented  in 
Volume  3,  Section  2.  The  exit  conditions  from  the  intermediate  zone 
start  the  calculations  for  this  zone  and  integration  of  the  nitric 
oxide  rr.te  equation  can  be  performed  to  the  end  of  the  combustor. 

In  practice,  the  state  of  uniform  mixing  as  idealized  above  is 
rarely  accomplished,  but  as  the  nitric  oxide  formation  rate  is  re¬ 
ducing  rapidly  in  this  zone  of  the  combustor,  it  was  anticipated  that 
the  assumption  would  have  little  significance  upon  the  levels  of  nitric 
oxide  predicted  to  exist  in  the  combustor  exhaust. 

3-3  MODEL  INPUT  REQUIREMENTS 

In  order  to  obtain  a  prediction  of  the  nitric  oxide  emissions 
level  from  the  model  described  above,  the  following  data  must  be  speci¬ 
fied: 


Kinetic  Constants:  The  reaction  scheme  for  nitric  oxide 

!  .  * 

formation  rate  described  in  Section  3*1  has  six  intermediate 
steps.  Each  step,  of  course,  has  its  own  characteristic 
values  for  A,  n,  and  E  and  these  must-be  specified  in  the 
input  data  (see  Table  2  for  values  of  these  parameters). 
Combustor  Operating  Conditions: 

T  ■  -  Inlet  temperature 

P  -  Operating  pressure 

-  Hass  mean  equivalence  ratio  in  the  primary  zone 
Combustor  Dimensions: 

Vp  -  Volume  of  the  primary  zone 

.  -  Length  of  the  intermediate  zone  (from  the  exit  of 
the  primary) 

Xr.,_  -  Distance  from  the  primary  zone  exit  to  the  exit  from 
tNP 

combustor  liner 

R  -  Radius  of  the  liner  at  position  X  for  all  positions 
downstream  of  the  primary  zone  exit  (for  annular 
combuscors  the  inner  and  outer  radii  must  be  specified) 
Mixing  and  Combustion  Parameters  for  the  Primary  Zone: 

Sq  -  Deg-ee  of  mixedness  i'n  the  prfmary  zone 

<8  -  Fraction  of  fuel  which  burns  in  primary  zone 

Air  Distribution: 

M  -  The  fraction  of  the  total  mass  flow  rate  of  air  mixed 
ax 

with  the  combustor  flow  stream  at  axial' posi tion  X 
(calculation  method  given  in  Volume;  3,  Appendix  XV:/ 
Mixing  Parameters  for  the  Intermediate  Zone: 

C^  -  A  parameter  specifying  the  rate  of  intermixing  between 
flow  elements  i..  the  intermediate  zone  (see  Volume  3, 
Section  2) 

Aj  -  Defines  the  rate  at  which -if-  collapses  in  the  inter¬ 
mediate  zone  according  to  the  relationship.  = 

O-o0  -  (x/xl)a1) 


A„,A  -  Define  the  .ate  at  which  unburned  fuel,  Mr  ,  is  burned 
2  3  •  fu 

in  the  intermediate  zone  by  M^^  =  >”fu  0  -  A^(X/X^)  2) 
Va'.ues  of  the  three  A  parameters  can  be  estimated  from 
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turbulent  mixing  phenomena  (see  Volume  3>  Section  2)  to  he 
given  by 

Aj  =  0-5 

a2  =  0.5 
=  1  *  ^ 

The  value  of  needs  to  be  determined,  but  as  will  be  shewn 
in  Section  4,  nitric  oxide  predictions  are  very  insensitive 
to  the  value  prescribed  and  a  value  of  5  is  taken  as  the 
nominal  value-  The  effects  of  changes  in  values  for  Aj ,  A^, 
and  are  also  determined  in  Section  4. 

and  Sq,  these  input  parameters 
can  be  determined  directly  for  the  case(s)  of  interest.  The  methods 
of  calculating  the  first  two  variables  are  described  in  Volume  3, 
Appendix  XVI-  The  value  of  Sq  must  be  determined  independently  for  each 
combustor  design  to  be  analyzed-  At  the  present  time,  Sq  cannot  be 
evaluated  directly  from  combustor  design  characteristics,  but  must  be 
estimated  indirectly  through  comparisons  cf  predicted  and  measured  NO 
emission  rates.  By  this  means,  the  value  of  has  been  found  to  lie 
within  the  range 

0-3  <  S  <.  0-7 
o 

for  combustors  of  conventional  design. 

3.4  SUMMARY  OF  MODEL  OUTPUT 

The  model  described  above  has  been  programmed  for  solution  by  digi¬ 
tal  computer  (see  Volume  3)  and  the  results  from  a  typical  test  case 
have  the  following  general  form: 

a.  In  the  primary  zone  the  flow  is  considered  as  a  series  of 
elements  each  bounded  by  discrete  steps  in  equivalence  ratio 
over  a  range  of  values  from  0.0  to  approximately  3-0-  The 
results  specify  the  mass  fraction  and  time-average  nitric 
oxide  concentration  for  each  element  at  the  primary  zone  exit 
together  with  the  total  mass  average  values  of  nitric  oxide, 
temperature,  and  density.  The  characteristic  residence  time 


Except  for  the  variables  (p  ,  M  , 
r  ^p  ax 
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and  mass  ave.rage  equilibrium  values  of  carbon  monoxide, 
carbon,  and  hydrocarbons  are  also  computed, 
b.  The  intermediate  zone  results  are  identical  to  those  com¬ 
puted  for  the  primary  zone  and  are  calculated  at  discrete 
steps  taken  downstream  of  the  primary  zone  exit.  The  com¬ 
putations  continue  until  the  rate  of  production  of  nitric 
oxide  becomes  insignificant  when  compared  to  the  earlier 
formation  rates.  If  this  occurs  in  the  intermediate  zone, 
one  final  computation  is  made  to  determine  the  conditions  at 
the  exit  of  the  combustor;  otherwise  the  computation  continues 
into  the  dilution  zone  where  only  mass  average  values  are 
computed  due  to  the  assumption  that  all  gases  are  well-mixed 
in  this  zone. 
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4.  PARAMETRIC  STUDY 

Those  parameters  which  describe  the  combustor  dimensions  and 
operating  conditions,  groups  b  and  c  in  Section  3*3,  are  reasonably 
well  defined  for  any  given  combustor,  but  the  groups,  a,  d,  e,  and 
.  f  contain  parameters  that  are  not  only  less-well  defined,  but  also 
are  essentially  internal  to  the  computing  procedure.  A  division 
was  made  between  them  therefore,  and  the  parametric  study  was  first 
concerned  with  the  latter  four  groups  in  what  is  termed  a  "modal 
sensitivity  analysis".  The  objectives  of  this  sensitivity  analysis 
were: 

a.  To  determine  how  significant  a  role  these  parameters  played 
ir.  predicting  NOx  emissions. 

b.  To  determine  the  best  values  to  use  for  these  parameters 
during  the  testing  of  the  group  b  and  c  variables. 

The  second  phase  of  the  analysis  is  termed  the  "design  sensi¬ 
tivity  analysis"  and  is  concerned  with  more  practical  aspects  such 
as  the  effect  of  the  size  of  the  primary  zone,  the  mean  equivalence 
ratio  in  the  primary  zone,  and  the  engine  pressure  ratio  upon 
the  levels  of  nitric  oxide  emission.  Predictions  obtained  from 
this  second  study  are  compared  with  experimental  data  published 
soon  after  the  study  was  conducted. 

4.1  MODEL  SENSITIVITY  ANALYSIS 

Ir.  order  to  make  full  use  of  the  results  obtained  from  the 
model  sensitivity  analysis,  the  selected  combustor  dimensions  and 
operating  conditions  corresponded  to  those  of  one  of  the  two  com¬ 
bustors  tested  experimentally  in  the  following  phase  of  this  study 
(see  Table  1). 

Nominal  values  were  chosen  for  the  mixing  parameters  and  each 
was  varied  systematically  to  at  least  one  value  higher  and  one  value 
lewer  than  nominal  while  all  other  variables  were  maintained  constant; 
thus,  Aj  nominal  was  set  equal  to  0.5,  but  calculations  were  also  made 
with  values  equal  to  0.2  and  1.0.  Values  assumed  for  the  other  mixing 
parameters  are  listed  in  Table  1,  together  with  those  values  used  to 
represent  the  combustor  and  its  operating  conditions. 
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Four  sets  of  data  were  selected  for  the  kinetic  rate  constants. 

The  sources  of  these  data  are  discussed  in  Volume  2,  Section  2,  and 
in  Reference  9*  Two  sets  of  data  were  arbitrarily  selected  as  those 
which  would  predict  the  lowest  and  highest  emissions  of  nitric  oxide. 
The  third  set  was  used  in  Reference  10  for  piston  engine  emission 
work.  The  fourth  set  was  selected  because  it  was  the  result  of  a 
comprehensive  search  of  the  literature  (Ref  11).  All  four  sets  of 
values  for  the  reaction  rate  data  are  given  in  Table  2.  Additional 
calculations  were  made  to  determine  the  importance  of  reactions  4  to 
6  (see  Section  3*1)  by  setting  their  rate  terms  all  equal  to  zero  for 
'each  of  the  four  sets  of  data. 

The  results  obtained  from  this  sensitivity  analysis  are  summarized 
graphically  in  Figures  2,  3»  and  4  and  Table  3,  are  discussed  below. 

4.1.1  Mixedness  of  the  Primary  Zone,  SQ 


Model  predictions  (see  Figure  2,  lower  figure)  show  that  for  the 
combustor  under  consideration,  the  value  assumed  for  the  mixing  condi¬ 
tions  in  the  primary  zone  does  not  play  a  very  significant  role  in  the 
level  of  NO  produced.  Also  shown  in  Figure  2  however,  are  other  results 
obtained  from  a  more  detailed  study  of  the  effect  of  the  mixedness." 

From  this  latter  data,  it  can  be  observed  that ^  =  0.8  represents 

a  low  sensitivity  operating  condition,  whereas  at  {$>£  -  1.0,  large 

reductions  in  NO  emissions  are  predicted  as  S  increases.  This  varia- 

o 

tion  in  sensitivity  may  be  explained  by  a  careful  consideration  of  the 
physical  significance  of  Sq  and  the  nature  of  the  NO  kinetics. 


These  results  were  produced  during  the  model  development  phase  of  the 
program  and  were  described  in  Reference  9-  The  combustor  geometry  and 
operating  conditions  used  were  somewhat  different  than  those  used  in 
the  parametric  analysis. 

$  g  represents  the  effective  mass  mean  equivalence  ratio  in  the  primary 
zone  when  account  is  taken  of  the  actual  fraction  of  fuel  burned  in  the 
zone  {0  )..  This  parameter  will  be  the  one  discussed  in  this  section, 
rather  than  ^  ,  the  mass  mean  equivalence  ratio  in  the  primary  zone  cal- 
qulated  on  the  assumption  that  all  fuel  is  burned  in  the  zone  ( 0  =  1), 
(p  g  relates  more  directly  to  the  nitric  oxide  formation  characteristics. 
The  two  are  simply  related  according  to  the  equation  0  _  =  <p  . 
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Mixing  in  the  primary  zone  is  never  complete  at  the  time  of  com¬ 
bustion  and  consequently,  elements  burn  with  different  equivalence 
ratios.  The  value  of  Sq  represents  this  "spread"  in  equivalence  ratios 
relative  to  the  mass  mean  value,  and  Sq  =  0  corresponds  to  the  perfectly 
mixed  case.  Clearly,  as  increases,  the  fraction  of  the  gases  in  the 
primary  zone  that  burn  at  the  mass  mean  equivalence  ratio  decreases. 

If  this  mass  mean  value  corresponds  to  that  at  which  the  nitric  oxide 
formation  rate  is  maximum,  then  the  total  nitric  oxide  formed  will 
decrease.  On  the  other  hand,  if  the  formation  rate  is  very  low  at  the 
mass  mean  value,  then  as  mixing  quality  decreases  (Sq  increases),  some 
elements  of  mass  will  burn  to  produce  a  higher  rate  of  formation  of 
nitric  oxide  and  the  emission  levels  will  increase. 

Computed  results  show  that  essentially  all  nitric  oxide  is  formed 
within  the  range  0.75<-  0  ^  <  1.25  for  the  combustor  inlet  conditions 
studied  (see  Table  1),  and  that  the  maximum  formation  rate  for  nitric 
oxide  corresponds  to  an  effective  equivalence  ratio  between  0.95  and 
1-0.  The  formation  rates  are  found  to  decrease  rapid’y  as  <f>  ^  decreases 
or  increases  from  the  0-95  value.  It  is  not  surprising,  therefore,  th3t 
the  variation  in  sensitivity  of  NO  emission  levels  to  changes  in  primary 
zone  mixing  characteristics  and  the  effective  mass  mean  equivalence 
ratio  are  as  shown  in  Figure  These  results  are  highly  significant, 
because  they  demonstrate  that  at  .mixedness  values  greater  than  0-3, 
the  NO  emissions  become  almost  independent  of £  for  the  range  of 
<|)  £  examined.  Thus,  large  reductions  in  emissions  will  only  occur  with 
changes  in  primary  zone  equivalence  ratio  if  the  primary  zone  conditions 
are  almost  homogeneous.  The  results  of  this  study  clearly  indicate 
the  futility  of  attempting  to  predict  nitric  oxide  emission  levels 
unless  some  account  be  taken  for  the  inhomogeneities  that  are  present 
in  all  gas  turbine  combustors. 

if.  1.2  Primary  Zone  Fuel-Burning  Fraction.  JL 

The  fraction  of  fuel  burned  in  the  primary  zone,/S  ,  was  computed 
from  the  primary  zone  combustion  efficiency  correlation  in  Figure  9- 
The  nominal  value  for  the  combustor  conditions  studied  is  0-8,  and 
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values  of  &  from  0.6  to  1.0  were  used  in  the  model  sensitivity  analy¬ 
sis.  As  shown  in  Figure  3a,  variations  in  /3  have  a  strong  effect  on 
NO  emissions,  and  a  reversal  in  the  sensitivity  occurs  in  the  range 
covered.  The  strong  effect  of  (3  i:.  attributable  to  its  effect  on  mean 
primary  zone  temperature  which,  in  turn,  has  a  strong  effect  on  the 
rates  of  NO  formation.  The  reversal  in  the  /3-  sensitivity  curve  was 
not  explored  in  detail,  but  is  suspected  to  result  from  increases  in 
temperature  and  NO  formation  in  the  intermediate  zone  as  £3  decreases. 

1  •  3  Air  Distribution  Characteristics 

The  air  distribution  characteristics  are  essentially  defined  by 
specifying  the  function, 

Ma  -  Ma(X) 

where  is  taken  to  be  not  the  total  air  within  the  liner,  but  rather, 
the  air  that  is  fully  mixed  with  the  main  combustor  flow  stream.  The 
method  used  to  calculate  these  distribution  characteristics  is  based 
upon  the  mixing  laws  discussed  in  Section  3.2.2  and  is  described  in 
Volume  3,  Appendix  XVI.  Variation  in  mixing  rates  may  have  signifi¬ 
cant  effects  upon  the  levels  of  nitric  oxide  predicted  by  tii„  odel, 
and  were  evaluated  as  follows. 

The  computer  results  showed  that,  for  the  combustor  conditions 
simulated,  the  nitric  oxide  formation  process  was  essentially  frozen 
in  less  than  one  radius  distance  downstream  of  the  primary  zone  exit. 
This  behavior  reflects,  of  co  rse,  the  drop  in  temperature  in  the  com¬ 
bustor  that  occurs  dec  to  the  audition  of  air  at  the  liner  wall.  Over 
the  one  radius  travel,  the  rate  of  air  addition,  dH  /dX,  was  nearly 

d 

constant,  so  the  effect  of  mixing  rate  upon  the  nitric  oxide  emission 
predictions  was  determined  by  selecting  two  other  higher  values  for 
dMa/dX  (the  upper  value  approaching  the  instant-mixing  assumption), 
and  predicting  new  values  for  the  emission  levels.  The  results,  shown 
in  Figure  3b,  indicate  that  the  nitric  oxide  level  predicted  can  vary 
by  approximately  20  per  cent  beiow  the  value  calculated  by  the  mixing 
law. 
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4-1.4  Mixing  Parameters  of  the  Intermediate  Zone 


Mixing  of  the  primary  zone  combustion  products,  the  unburned  fuel, 
and  the  additional  air  entering  at  the  liner  walls  is  assumed  in  the 
model  to  follow  turbment  mixing  laws  and  to  be  approximately  propor¬ 
tional  tf'  the  square  root  of  the  fractional  distance  traveled  down  t^e 
intermediate  zone.  The  sensitivity  of  the  model  to  this  assumption 
was  checked  by  varying  the  exponent  to  values  above  and  below  the  value 
of  one-haif  and  the  results  are  considered  below. 

The  exponent  represents  the  mixing  characteristic  for  the  com¬ 
bustion  products  and  the  results  given  in  Figure  4  show  the  effect  of 
changing  its  value  about  the  nominal  value  of  one-half.  Nitric  oxide 
levels  are  shown  to  increase  from  -10  to  +10  per  cent  about  this 
nominal  value  as  A^  increases  from  0.2  to  a  high  value  of  1.0.  An 
increase  in  the  value  of  Aj  corresponds  to  a  decrease  in  the  mixing 
rate  and,  thus,  a  decrease  in  the  rate  of  reduction  of  mass  in  the  flow 
elements  with  near-stoichiometric  mixture  jatios.  Hence,  a  reduced 
mixing  rate  is  accompanied  by  ar  increase  in  total  NO  formation  for 
a  system  where  ^  <  1.0.  The  20  per  cent  variation  in  NO  actually 
represents  a  muc.n  larger  fraction  of  the  NO  formed  in  the  intermediate 


zone. 

The  exponent  A2  and  coefficient  A^  together  define  the  rate  of 
burn-Up  of  the  fuel  which  leaves  the  primary  zone  unburned  (see  Fig  4). 
These  parameters  effectively  determine  the  average  axial  temperature 
profile  in  the  early  stages  of  the  intermediate  zone  in  a  manner  which 
is  also  related  to  the  local  rate  of  addition  of  air  at  the  liner  walls 
For  the  nominal  case.  $  -  =  fi-8  and  0  =0.8,  an  increase  in  the  rate 

of  fuel  burning  (A^  smaller,  A^  larger)  results  in  the  maintaining  of 
a  high  temperature  condition  in  the  early  stages  of  the  zone.  As  s'sown 
in  Figure  4,  however,  the  NO  emission  rate  for  this  case  does  not  vary 
greatly —  only  from  -10  to  +20  per  cent  of  the  nominal  value.  However, 
as  with  the  effect  of  Aj ,  this  30  per  cent  variation  of  the  total  NO 
emission  rate  represents  a  large  fraction  of  the  NO  formed  in  the  inter 
mediate  zone. 

The  final  mixing  parameter  ,  introduced  for  this  zone  defines 
the  r3te  of  intermixing  between  the  elements  of  various  equivalence 
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ratios  and  nitric  oxide  concentrations-  It  specifies  the  amount  of 
mass  that  leaves  an  element  due  to  intermixing  during  an  incremental 
step  down  the..zone  to  the  total  mass  flow  rate  through  the  element 
(see  Volume  3,  Section  2).  The  greater  the  value  of  C^,  the  more 
vigorous  the  mixing  between  the  elements-  Its  main  effect  is  to 
modify  the  nitric  oxide  concentration  and  so  influence  the  local 
formation  rate.  The  net  result  of  all  mixing  is  to  change  element 
concentrations  towards  the  average  value,  thus  certain  elements  will 
increase  in  concentration  and  others  reduce  in  concentration.  An 
increase  in  concentration  will  reduce  the  nitric  oxide  formation  rate 
and  vice  versa.  It  is  not  surprising  therefore,  to  find  that  the  model 
predicts  that  nitric  oxide  emissions  are  also  quite  insensitive  to  the 
value  given  this  parameter  (see  Fig  k) ,  as  the  effects  are  self-com¬ 
pensating. 

4.1.5  Kinetic  Rate  Data 

The  results  obtained  by  using  the  alternative  sources  of  reaction 
rate  data  described  previously,  are  summarized  in  Table  3-  Nitric 
oxide  emission  levels  can  be  observed  to  vary  by  a  factor  greater  than 
two,  depending  upon  the  set  of  data  used.  This  "parameter"  is  shown 
to  exert  the  greatest  influence  upon  the  model  predictions.  The 
results  also  show  that  reactions  4,  5,  and  6  exert  a  negligible 
influence  upon  the  results  for  all  cases  tested. 

4. 1 . 6  Conclusions 

As  stated  earlier,  the  purpose  of  the  model  sensitivity  analysis 
was  to  evaluate  the  sensitivity  of  the  model  predictions  to  variations 
in  model  parameters  (as  distinguished  from  design  parameters).  With 
this  purpose  in  mind,  the  following  conclusions  can  be  drawn  from  the 
model  sensitivity  analysis  for  the  case  tested: 

a.  The  degree  of  mixedness,  S  ,  and  the  fuel-burned  fraction, /3  , 

o 

assigned  to  the  primary  zone  have  significant  effects  upon  the 
predicted  nitric  oxide  emissions. 
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b.  Mixing  rates  in  the  intermediate  zone  affect  the  amount  of 
NO  formed  in  this  zone  to  a  substantial  degree-  This  effect 
is  most  pronounced  with  the  rates  of  mixing  of  elements  of 
differing  fuel  fraction;  and  mixing  of  unburned  fuel  and 
secondary  air.  For  the  combustor  design  and  operating  condi¬ 
tions  investigated,  most  NO  is  formed  in  the  primary  zone. 
Therefone,  intermediate  zone  mixing  effects  are  considered  to 
be  of  secondary  importance  with  respect  to  total  NO  emission 
rate. 

c.  The  predicted  NO  emission  rate  is  very  sensitive  to  the 
reaction  rate  constants  used  in  the  NO  formation  model.  It 
is  shown  that,  depending  on  the  data  selected,  predicted 
nitric  oxide  levels  can  vary  by  a  factor  of  greater  than 
two.  Some  variation  in  predicted  NO  emissions  was  expected 
since  a  wide  selection  of  rate  constants  was  used.  However, 
as  indicated  in  Table  3,  the  greatest  discrepancy  occurs  be¬ 
tween  the  twc  sets  of  constants  which  were  considered  to  be 
most  reliable  for  this  application  (l CODE'S  1  and  2). 

4.2  DESIGN  SENSITIVITY  ANALYSES 

Conditions  simulated  in  the  design  sensitivity  analysis  were 
those  experienced  by  a  combustor  when  operating  under  conditions 
corresponding  to  the  idle,  take-off,  cruise,  and  approach  modes  of  a 
mode.n  airliner.  The  three  variables  considered  to  be  the  most  signi¬ 
ficant  and  studied  in  detail  were  as  follows: 

a.  The  mean  residence  time  :n  the  primary  zone,  ^  - 

b.  The  mean  equivalence  ratio  in  the  primary  zone,  . 

c.  The  engine  pressure  ratio,  EPR. 

Throughout  this  analysis,  the  condition  was  imposed  that  the 
thrust  level  be  constant  for  each  operating  mode  despite  changes  in 
these  p-ime  variables. 


4.2.1  Mean  Primary  Zone  Residence  Time 

if  all  conditions  are  constant,  then  the  primary  zone  residence 
time  is  directly  proportional  to  the  volume  of  the  primary  zone,  V  . 


The  effect  of  residence  time  upon  NO  emission  characteristics  was 
determined  at  each  of  the  four  operating  modes  simply  by  assuming 
that : 

V  -  0.7  x  Nominal  Value 
P 

V  *=  1.3  x  Nominal  Vaiie 
P 

The  results  are  given  in  Figure  5*  Ail  show  a  strikingly  similar 
variation  in  nitric  oxide  emission  levels  with  the  level  increasing 
by  approximately  1.2  lbs  N0/1Q00  lbs  of  fuel  for  all  modes  except 
the  idle  mode. 

It  is  interesting  to  observe  that  the  nitric  oxiJe  emissions 
do  not  increase  in  direct  proportion  to  the  primary-i-nyie  residence 
time  particularly  at  the  operating  conditions  which  correspond  to  the 
high  levels  of  emission.  This  result  can  be  attributed  to  the  fact 
that,-  in  the  combustor  studied,  a  significant  fraction  of  the  total 
NO  formed  is  found  in  the  intermediate  zone.  In  the  analysis,  the 
intermediate  zone  conditions  were  not  changed  so  that  the  NO  fraction 
formed  there  also  did  not  change.  The  NO  fraction  formed  in  the  primary 
zone  varies  almost  in  direct  proportion  to  primary  zone  residence  time, 
^  •  Thus,  the  total  NO  emission  level  does  not  vary  in  direct  pro¬ 
portion  to  H  . 

4.2.2  Mean  Primary  Zone  Equivalence  Ratio 

A  change  in  the  primary  zone  equivalence  ratio  will  charge  the 
residence  time  if  the  primary  zone  volume  is  kept  constant  and  will 
also  produce  a  change  in  the  over-all  equivalence  ratio  (hence,  thrust 
level)  unless  the  airflow  characteristics  are  modified.  For  che  oper¬ 
ating  conditions  simulated  in  this  study,  residence  time  was  specified 
as  constant  for  each  mode  and  equal  to  the  nominal  value.  The  constant 
thrust  condition  imposed  at  the  onset  of  the  calculations  was  assumed 
to  be  satisfied  if  the  total  mass  flow  rates  of  the  fuel  and  air  were 
maintained  constant.  For  richer  primary  zones,  therefore,  less  air 
was  assumed  to  enter  the  primary  zone,  and  the  air  distribution  char¬ 
acteristics  were  modified  such  that  the  total  airflow  at  the  exit 
remained  constant  for  each  mode.  The  following  simple  relationship  was 
used: 
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M  =  M 
a  a 


i  -  (.  -  i/a)  <VA 


nominal  L 

where  K  =  total  airflow  rate 
a 

R  =  the  ratio  of  the  tested  equivalence  ratio  to  the  nominal 
value 

X  =  the  axial  position  corresponding  to  the  value  of 
=  the  end  of  the  combustor. 

This  method  was  used  to  determine  the  effect  that  the  primary  zone 
equivalence  ratio  has  upon  NO  emission  levels  at  the  values  of  the 
ratio  other  than  the  nominal  one.  These  corresponded  to 

3  =  (  <£  )  X  0.8 

”  p  ”  p  nom 

x  =  (  3  )  x  1.2 
p  t  p'  nom 

and  all  four  operating  modes  wers  tested.  The  results  are  shown  in 
Figure  6. 

As  could  be  anticipated,  in  those  cases  in  which  the  change  results 

in  <J>  values  away  from  the  stoichiometric  condition,  NO  emissions 

reduce,  whereas  when  ©  moves  toward  this  condition,  the  NO  emissions 

v  p 

increase.  This  effect  is  most  significant  for  the  cruise  and  approach 
operating  modes,  where  for  the  test  cases,  NO  emissions  are  shown  to 
increase  by  a  factor  of  three  over  the  range  of  primary  zone  equiva¬ 
lence  ratios  considered.  Special  note  must  be  made,  however,  of  the 
influence  of  the  assumed  value  of  the  primary  zone  mixedness  parameter, 
Sq.  The  nominal  value  used  in  these  calculations  corresponded  to  a 
value  of  0.2  (Table  1).  The  sensitivity  of  these  results  to  the  value 

of  S  as  shown  in  Figure  2. 
o  3 


^.2.3  Pressure  Ratio 

For  the  nominal  test  case,  the  engine  pressure  ratio  was  equal  to 
16.7:1  at  the  take-off  mode.  Operating  conditions  were  simulated  for 
the 

EPR  (take-off)  =  12:1  and 
EPR  (take-off)  =  25:1 

and  their  nitric  oxide  emissions  were  predicted  at  all  four  modes  of 
aircraft  operation. 
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Pressure  levels  within  the  combustor  aL  off-design  conditions 
were  assigned  to  be  proportional  to  the  EPR  and  were  computed  based 
upon  the  data  for  the  nominal  case.  Ccmbustor  inlet  temperatures 
were  determined  from  data  supplied  in  Hodge  (Fig  104,  Ref  i2). 

In  order  to  maintain  operating  thrust,  primary  zone  mean  resi¬ 
dence  time  and  equivalence  ratio  all  constant,  and  at  the  same  time 
vary  the  EPR,  new  combustor  dimensions  had  to  be  defined  for  the 
cases  of  EPR  =  12:1  and  25:1-  The  method  used  assumed  the  following: 

a.  That  the  combustors  operated  with  the  same  reference  velo¬ 
city,  vr,  as  for  the  nominal  ccmbustor,  where: 

vr  =  M//°  A  and 
M  total  airfiow  rate 
P  =  air  density  at  inlet 
A  =  cross-sectional  area  of  liner 

If,  as  an  approximation,  the  density  is  assumed  to  be  pro¬ 
portional  to  pressure,  P,  and  if  the  v.otal  mass  flow  is  fixed, 
then  this  criterion  implies  that  the  product  PA  is  constant. 
The  value  of  the  product  was  known  for  the  nominal  combustor, 
so  the  values  of  the  liner  area  could  be  computed  for  the 
two  combustors  operating  at  other  pressure  levels. 

b.  That  the  volumes  of  the  primary  zone  of  the  new  combustors 
were  such  that  the  residence  times  were  equal  to  that  of 
the  nominal  combustor.  This,  of  course,  is  an  artificial 
constraint  upon  the  combustor  dimensions,  but  it  was  applied 
so  that  the  effect  of  pressure  ratio  upon  nific  oxide  emis¬ 
sions  could  be  evaluated  in  isolation  from  the  other  prime 
variables.  Checks  were  made  to  ensure  that  the  combustors 
sized  by  this  approach  would  have  an  adequate  combustion 
efficiency  by  methods  detailed  in  Reference  4.  The  methods 
showed  that  for  the  EPR  values  tested,  ali  three  combustors 
would  perform  with  the  same  efficiency. 

The  effect  of  EPR  upon  nitric  oxide  emissions  rests  not  so  much 
upon  the  absolute  pressure  level  that  results,  as  upon  the  fact  that 
the  inlet  tempercture  to  the  combustor  and  flane  temperature  increase 
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with  the  EPR.  At  high  inlet  temperatures  of,  say,  greater  than 
800  deg  K,  changes  of  the  order  of  50  deg  K  have  significant  effects 
upon  the  predicted  nitric  oxide  emission  level.  The  results  given 
in  Figure  7  show  that  the  engine  compressor  pressure  ratio  has  a  far 
greater  effect  upon  emission  level  of  nitric  oxide  than  all  other 
variables  studied;  that  is  to  say  that  the  inlet  temperature  to  the 
combustor  tends  to  exert  the  greater  influence  upon  emission  charac¬ 
teristics.  This  applies  for  all  operating  modes  (except  idle),  and 
it  is  predicted  that  as  the  EPR  progresses  from  12  to  16-7  (the  nominal 
value)  and  to  25,  er  'ssion  levels  at  all  modes  increase  approximately 
in  the  ratios  1:2:^. 5*  In  practice,  this  dependence  is  not  so  large, 
as  other  factors  must  be  considered  when  the  value  of  EPR  is  changed 
and  a  discussion  of  these  factors  follows. 

4. 2.*»  Pressure  Ratio  Effects.  Specific  Fuel 

Consumption.  Specific  Thrust  and  Emissions 

A  change  in  EPR  will  result  in  a  change  in  both  the  thermo¬ 
dynamic  efficiency  of  the  engine  and  a  change  in  the  specific  thrust, 
i.e-,  the  lbs  thrust  per  lb  of  air  per  sec  that  flows  through  the 
engine-  In  general,  the  specific  thrust  is  relatively  insensitive 
to  EPR  at  the  conditions  tested  above  (i.e.,  a  fixed  exit  temperature 
from  the  combustor),  and  upon  this  observation  is  based  the  assumption 
stated  earlier  that  constant  thrust  corresponds  to  a  fixed  mass  flow 
rate.  But  clearly,  any  improvement  in  the  thermodynamic  efficiency 
must  result  in  less  fuel  being  burned  per  unit  time  to  achieve  the 
same  thrust  level,  and  hence  less  nitric  oxide  being  generated.  Some 
credit  should  be  allocated  to  the  higher  EPR  engines  for  this  fact. 

For  the  pressure  range  and  operating  conditions  tested  above,  however, 
this  credit  is  within  j0.05  of  the  nominal  value  and  can  therefore 
be  ignored. 

The  considerable  improvement  in  engine  performance  that  can  be 
obtain*  '  by  resorting  to  higher  EPR  values  is  only  realized  if  the 
maximum  cycle  temperature,  T  (that  is,  the  mean  temperature  at  the 
exit  of  the  combustor),  is  also  increased.  Such  a  change  does  further 
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improve  specific 

fuej  consumption,  but  has  a  much 

greater 

upon  the  specific 

thrust  as  is  shown  below: 

: 

? 

Specific 

Over -A 1 1 

.  Deg  K 

Thrust  Ratio: 

(Approx) 

1273 

1.0 

0. 17 

1473 

1 . 25  ! 

0.27 

1673 

1.45 

0.37  : 

1873- 

1.60 

0.47 

where  specific  thrust  ratioi  is  the  thrust  at  T  divided  by  that  at 
_  m  ; 

T  equal  to  1273  deg  K,  which  corresponds  approximately  with  the 
nominal  operating  point.  These  results  indicate  that  the  same  thrust 

i 

level  can  be  delivered  for* a  lower  mass  flow-  rate  of  air.  It  is 
•interesting  to  note  also  that  the  required  over-all  fuel-to-air  ratio 
never  exceeds  one  half  of  the  stoichiometric  requirement.  As  all  nitric 

i  J  • 

oxide  formation  reactions  are  essentially  frozen  at  ratios  less  than 
six  tenths  the  stcichiometr ic  value,  this  implies  that,  to  a  first  ' 
approximation,  the  nitric  oxide  concentrations  at  the  combustor  exit 
need  not  be  changed  at  all  by  the  increase  in  the  maximum  cycle  tem¬ 
perature.  Hence  some  credit  is  due  to  the  combustor  with  the  higher 
specific  thrust  as  the  total  nitric  oxide  emitted  per  unit  time  is 
decreased.  The  maximum  cycle  temperature  that  is  currently  in  common 
practice  corresponds  to  a  specific  thrust  ratio  of  about  1.5.  If  this 
correction  is  -applied  to  the  emission  ratios  given  at  the  end  of  the 
preyious  section,  then  the  predicted  nitric  oxide  emission  rates  are 
approximately  increased  by  1  =; 2 : 3  as  the  EPR  progresses  from  12  to 
16.7  to  25-  1 


A. 2. 5  Experimental  Data 

In  Figure  8,  measured  NOx  emission  data  are  plotted  in  the  same 

i  ■  ;  , 

manner  as  the  theoretical  predictions  'in  Figure  7-  The  data  have  been 
taken  from  Reference  13,  and  each  point  in  the  figure  is  an  average 
measurement  for  several  different  engines  of  the  same  model  or  family. 
A  comparison  of  Figures  7  and  8  indicates  that  the  effects  of  engine 
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pressure  ratio  on  NO  emission  rate  as  predicted  by  the  model  are 
substantially  confirmed  by  the  measured  data.  The  main  discrepancy 
occurs  at  the  idJe  condition  where  the  model  predicts  essentially 
zero  emissions. 


^.2.6  Conclusions 


The  following  conclusions  can  be  made  from  the  results  obtained 
in  the  design  sensitivity  analysis: 

a.  The  model  does  predict  nitric  oxide  emissions  which  are  in 
accord  with  measured  values. 

b.  The  dominating  influence  upon  emissions  level  is  the  engine 
pressure  ratio,  but  large  reductions  in  nitric  oxide  emis¬ 
sions  appear  to  be  feasible  for  any  given  operating  pressure 
by  modifying  the  primary  zone  operating  conditions. 


5.  THE  EXPERIMENTAL  PROGRAM 
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5.1  TEST  OBJECTIVES 

The  results  described  in  the  previous  sectiors  indicate  that, 
although  the  present  computer  mod-»l  does  predict  nitric  oxide  emission 
levels  that  are  in  reasonable  accord  with  experimental  values,  there 
is  some  question  concerning  the  accuracy  in  calculation,  or  assumed 
value,  of  some  of  the  parameters  used  in  the  computation.  The  main 
objective  of  the  experimental  test  program  was  to  determine  the  nitric 
oxide  emission  levels  from  some  gas  turbine  combustor?  under  closely 
controlled  test  conditions  so  that  the  relevance  of  these  questions 
could  be  resolved  and  confidence  in  the  model  increased. 

At  the  same  time  it  was  planned  to  collect  emission  levels  of 
other  pollutant  species  as  they  varied  with  the  test  operating  conditions 
and  so  define  the  emission  characteristics  for  these  species.  This 
data  could  then  be  U6ed  to  indicate  the  effect  that  nitric-oxide 
reducing  changes  in  combustor  design  might  have  upon  the  levels  of 
thase  emissions. 

5-2  TEST  PLAN 

It  is  concluded  from  the  foregoing  analysis  that  for  a  fixed  combustor 
and  operating  condition  the  following  parameters  are  the  most  significant 
in  terms  of  their  effect  upon  nitric  oxide  emission  levels: 

a.  The  fuel-burning  fraction  in  the  primary  zone,fi. 

effect  of  imperfect  m  xing  in  the  primary  zone,  S  . 

o 

c.  The  air  distribution  characteristics,  dM  /dx. 

d 

d.  The  kinetic  rate  constant:. 

Of  these  four  parameters,  only  three  are  variables  as  there  can 
be  only  one  accurate  set  of  kinetic  data.  Due  to  the  o-esent  uncertairtv 
in  the  relative  accuracy  of  the  cat  3  available  in  the  literature, 
however,  the  correct  one  must  be  se  ected  based  upon  the  experimental 
results.  Clearly,  the  test  plan  mu>t  aim  to  produce  conditions  for 
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which  significant  variations  in  the  three  variables  occur.  The  test 
plan  adopted,  therefore,  was  to  test  two  combustors  developed  by 
different  design  criteria  (i.e. ,  different  manufacturers)  and  the 
following  tests  were  made. 


5.2.1  Test  Series  1  -  Combustor  A 

Exhaust  concentrations  only  were  measured  in  the  first  test  series. 
Concentrations  of  nitric  oxide,  carbon  monoxide,  and  unburned  hydro¬ 
carbons  were  recorded  as  the  operating  conditions  of  the  combustor  were 
systematically  varied  to  cover  the  following  ranges  of  operating  condi¬ 
tions;  where  P  =  inlet  pressure  in  psia,  T  =  inlet  temperature  in  deg  K, 
APR  =  air-fuel  ratio. 


p 

1 

AFR 

Comments 

85 

600 

70 

Residence  time 

85 

700 

50  to  140 

varied  for 

85 

800 

60 

same  test  conditions 

115 

600 

55  to  140 

Residence  time 

115 

700 

50  to  120 

at  design 

115 

850 

50  to  90 

values 

40 

400 

40  to  80 

Residence  time 

at  design  value 

est  Series 

2  -  Combustor  B 

the  second 

series. 

measurements  were 

taken  both  at  the  exit 

:  the  combustor  and 

at  axial  stations 

•within  the  combustor  a 

owing  conditions: 

P 

T 

AFR 

Comments 

85 

600 

60  to  130 

Exhaust  concentrations 

85 

700 

55  to  130 

only,  residence  time 

85 

800 

65  to  115 

at  design  value 

70 

600 

88 

Internal  probing 

Measurements  were  taken,  within  the  combustor  of  radial 
profiles  of  nitric  oxide,  carbon  monoxide,  temperature,  and 
pressure  at  two  axial  stations  located  approximately  0.5 
diameters  and  1.0  diameters  downstream  of  the  primary 
zone  exit. 
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5-3  EXPERIMENTAL  EQUIPMENT 

5*3-1  The  Combustor  Testing  Facility 

The  experimental  test  arrangement  employed  is  shown  schematically 
in  Figure  10  and  photographically  in  Figure  11.  It  comprised  an  air 
supply,  an  air  preheater  in  the  form  of  a  pebble-bed  heat  exchanger, 
an  air  mixer,  a  water-cooled  tailpipe,  and  a  high  temperature  regula¬ 
ting  valve.  Control  was  effected  remotely  in  a  test  rocm  located 
about  15  ft  from  the  combustor  (see  Fig  12).  Airflow  rates  for  most 
tests  were  within  the  range  of  4  to  7  lbs  per  second  and  combustor 
inlet  temperature  was  adjusted  by  diverting  different  percentages  of 
air  through  the  pebble-bed.  Temperature  traverse  quality  at  the  exit 
of  the  mixer  was  within  +15  deg  K  of  the  nominal  value,  and  the  com¬ 
bustor  entrance  was  located  just  downstream  of  the  mixer  exit.  Casings 
were  manufactured  co  the  dimensions  experienced  by  the  combustors 
during  in-flight  conditions  to  ensure  that  the  air  distribution  charac¬ 
teristics  in  operation  during  testing  matched  design  conditions. 

5-3-2  Sampling  Methods 

Two  sampling  methods  were  employed.  The  first  system,  shown 
in  detail  in  Figure  1 3 »  collected  exhaust  gases  from  the  exit  plane 
of  the  combustor  through  a  diametral  probe.  The  probe  contained  nine 
holes  of  0.040  inches  in  diameter  spaced  to  give  equal  area  sampling. 

The  probe  itself  was  water-cooled,  but  the  20  ft  of  sample  line  from 
the  probe  to  the  instrumentation  room  was  heated  to  about  470  deg  K 
to  prevent  condensation  of  the  hydrocarbons.  After  hydrocarbon  sampling, 
the  gases  passed  through  a  series  of  coolers,  traps,  and  filters  to 
remove  the  water  and  the  carbon. 

The  other  probe  system  employed  was  used  for  probing  the  interior 
of  the  second  combustor  tested.  A  stainless-steel  probe  was  water- 
cooled  with  a  high  pressure  water  supply  system  and  entered  the  combustor 
normal  to  the  direction  of  the  gas  flow.  The  probe  lip  consisted  of  a 
short  length  of  right-angled  silica  tube  which  was  inserted  well  into 
the  metal  body  of  the  probe.  Silica  was  used  to  limit  the  possibility 
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of  obtaining  erroneous  readings  for  nitric  oxide  concentration  levels 
if  the  sampling  region  probed  were  to  be  fuel-rich  (see  Ref  14).  The 
exit  of  the  probe  was  then  connected  to  the  sample-conditioning  train 
used  for  the  exhaust  probe  (see  Fig  13)- 

5-3*3  Instrumentation 

Measurements  of  airflow  rate,  pressure,  temperature,  and  velocity 
were  all  made  using  conventional  instrumentation  techniques.  Pollutant 
specie  concentrations  were  measured  as  follows: 

a.  Nitric  Oxide.  Four  techniques  were  used  for  measuring  nitric 
oxide  concentrations  but  only  one  was  used  consistently  through¬ 
out  the  course  of  the  test  program.  This  was  a  nondispersive 
infrared  method  which  employed  a  Grubb  Parsons  SB2  instrument 
with  a  20  cm  cell  and  a  specially-developed,  integral,  filter 
system  to  minimize  the  response  to  water  and  hydrocarbons.  It 
was  found,  however,  that  the  measurements  taken  using  this 
method  were  invalidated  due  to  any  water  in  the  sampled  gas. 

The  problem  was  eliminated  by  use  of  an  additional  drying  tube 
packed  with  calcium  sulphate.  The  system  sensitivity  to  any 
unburned  hydrocarbons  in  the  sample  gas  was  always  suspect 
until  it  was  checked  against  a  Thermo  Electron  chemiluminescent 
analyzer.  This  analyzer  was  employed  continuously  over  the 
last  one-third  of  the  test  program  and  showed  that  excellent 
agreement  could  be  obtained  between  the  two  techniques  for 
all  nitric  oxide  levels  greater  than  about  15  ppmV  nitric 
oxide.  These  emission  levels  corresponded  to  simulated  idle 
conditions  in  the  engine  which  also  had  very  high  levels  of 
unburned  hydrocarbons  (greater  than  100  ppmV).  Even  under 
these  conditions,  the  chemiluminescent  analyzer  gave  readings 
which  were  no  lower  than  10  ppmV  below  the  ND1R  readings-  In 
operation,  both  methods  were  calibrated  with  pre-cal ibrated 
gases  of  200  and  100  ppmV  NO  in  nitrogen  and  linearity  was 
checked  by  controlled  dilution  of  these  gases.  Linearity  was 
found  to  be  better  than  required  in  both  cases  (-*-2  ppmV  NO). 
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Measurements  were  also  using  the  standard  Saltzman 
and  Phenol  Disulphonic  Acid  techniques  but  these  were 
restricted  mostly  to  test  conditions  corresponding  to  the 
idle  mode  and  before  the  cnemi lumi nescent  instrument 
became  available  for  use. 

fc.  Carbon  Monoxide.  Carbon  monoxide  was  also  measured  using 
the  NDIR  method  by  means  of  a  Grubb  Parsons  SB I  analyzer. 
This  analyzer  had  a  30  cm  cell  with  a  carbon  dioxide  filter 
attached. 

c.  Unburned  Hydrocarbons.  A  Perkin-Elmer  FI1  flame  ionization 
gas  chromatograph  was  used  to  measure  unburned  hydrocarbons 
and  the  instrument  was  calibrated  with  hexane. 

S.h  SUMMARY  OF  RESULTS 

A  summary  of  all  test  resul's  obtained  is  given  in  Tables  h 
and  5  and  these  will  be  presented  in  graphical  form  in  the  next 
section  where  they  are  compared  with  the  model  predictions. 
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6.  COMPARISON  BETWEEN  MODEL  PREDICTIONS  AND  EXPERIMENTAL  RESULTS 

6.1  MODEL  INPUT  DATA 

The  required  input  data  for  the  model  were  first 
determined  for  each  of  the  two  combustors  tested  experimentally 
and  according  to  the  f'inciples  described  in  Sections  3  and  4. 

6.1.1  Combustor  Characteristics 

The  greater  part  of  this  calculation  involved  the  determina¬ 
tion  of  the  air  distribution  characteristics  throughout  each  com¬ 
bustor  and,  in  both  cases,  the  distribution  characteristics  were 
based  upon  measured  data  supplied  by  the  manufacturers.  These 
data  specified  the  percentage  of  thr  total  air  supplied  to  the 
combustor  that  entered  at  each  hole  in  the  combustor  liner,  and  the 
objective  of  the  calculations  was  to  take  account  of  the  finite 
distance  lag  that  must  exist  before  the  entering  air  gets  entrained 
by  the  combustion  products  formed  in  the  upstream  sections  of  the 
combustor.  The  method  used  is  detailed  in  Volume  3,  Appendix  XVI, 
and  the  results  are  shown  graphically  in  Figure  14.  The  combustor 
dimensions  are  given  in  the  table  below,  and,  like  the  air  distribu¬ 
tion  characteristics,  the  values  were  considered  invariant  with 
combustor  operating  conditions: 

Combustor  Geometry 


Volume  of  Primary  Zone  (V  )  in^ 

Length  of  Intermediate  Zone  (X^)  in 
Length  of  Intermediate  and  Dilution  Zone  in 

Combustor  Radius  (R)  in 

6.1.2  Combustor  Operating  Conditions 

Five  variables  are  sufficient  to  describe  the  combustor  oper¬ 


Combustcr 


A 

8 

38.0 

88-5 

7-0 

7-0 

14.0 

12-5 

2.63 

3-25 

ating  conditions: 
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a.  T  -  inlet  temperature 

b.  P  Operating  pressure 

c.  p  -  Mass  mean  equivalence  ratio  in  the  primary  zone  (for =  1) 

d.  &  -  Fraction  of  fuel  entering  the  primary  zone  which 

r 

burns  in  the  zone 

e.  M  The  total  mass  airflow  rate  into  the  combustor. 

ax 

Conditions  a,  b,  c,  and  e  were  made  to  match  the  conditions 
tested  experimentally  and  the  value  of was  calculated  from  the 
correlation  given  in  Reference  15  and  also  shown  for  convenience 
in  Figure  9*  The  value  of  the  primary  zone  equivalence  ratio  is 
simply  related  to  the  over-all  air-fuel  ratio  in  a  manner  dictated 
by  the  fraction  of  air  that  enters  the  primary  zone,  and  for 
the  two  combustors  tested,  this  relationship  is  shown  in  Fig¬ 
ure  15*  Also  shown  in  this  figure  is  the  effective  primary 
zone  equivalence  ratio  chat  exists  when  account  is  taken  of  the 
primary  zone  fuel-burning  fraction,^,  for  one  particular 
operating  pressure.  It  can  be  observed  from  this  figure  that  the 
operating  conditions  in  the  primary  zones  vary  significantly  for  the 
two  combustors  tested  experimentally.  Combustor  A  is  calculated 
always  to  operate  3t  conditions  below  the  stoichiometric  condition, 
whereas  Combustor  B  is  shown  to  operate  fuel -rich  for  air-fuel  ratios 
less  than  65- 


6.1.3  Primary  Zone  Mixing  Parameter 

This  parameter,  Sq,  describes  the  degree  of  mixedness  in  the 
primary  zone  (see  Section  3*2.1).  Practically,  this  parameter 
can  have  a  range  of  values  from  perhaps  near  zero,  corresponding 
to  near  perfect  mixing,  to  about  one,  very  poor  mixing,  and  for  the 
purpose  of  comparing  the  model  predictions  with  the  experimental 
results,  it  was  considered  to  be  a  matching  parameter.  That  is  to 
say,  its  prescribed  value  was  varied  parametrically  and  the  effects 
of  so  doing  observed  in  the  predicted  results  and  compared  with  the 
experimental  results. 
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6.1.4  Kinetic  Constants 

The  reaction  scheme  for  nitric  oxide  formation  rate  described 
earlier  has  six  intermediate  steps,  each  with  the  general  form, 

k  -  ATn  exp(-E/RT) 

where  k  is  the  reaction  rate  constar-t.  Each  step,  of  course,  has 
its  own  characteristic  values  for  A,  n,  and  E  and  these  must  be 
specified  in  the  input  data.  Four  sets  of  such  data  are  presented 
in  Table  2  and  are  referenced  as  I  CODE  =1  to  4. 

At  the  start  of  this  study,  the  first  set  of  data  was  thought 
to  be  the  most  useful  as  it  had  been  successfully  used  to  correlate 
piston  engine  nitric  oxide  emissions  (Ref  10).  The  second  set  of 
data  (ICODE  =  2)  was  considered  the  most  accurate  as  it  is  the 
product  of  an  extensive  study  of  all  the  available  nitric  oxide 
kinetic  data  (Ref  11).  Results  obtained  with  the  latter  data 
predicted  higher  values  of  nitric  oxide  emissions  than  the  former 
and  this  set  of  data  has  been  used  to  obtain  the  predicted  results 
presented  in-this  section.  Results  obtained  with  the  first  sat  of 
data  were  found  not  to  matcn  the  measured  results  as  well  as  with 
the  second  set. 

Measured  First  NO  Predictions^ 

Value  Kinetic  Data  (ICODE) 

1  2 

ppmV  ppmV  ppmV 

132  70  145 

Some  typical  results  are  shown  in  the  above  table.  It  is  recommended 
that  the  data  referenced  by  the  code  ICODE  =  2  be  used  with  the  model 
presented  in  this  report. 

See  Table  2.  ICODE  =  1  corresponds  to  data  used  to  match  results  for 
piston  engine  and  (CODE  **  2  corresponds  to  the  data  selected  in  a 
comprehensive  survey  of  the  data  as  being  the  most  reliable  (see  Ref  11). 
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6.2  MODEL  VERIFICATION 

The  two  combustors  employed  in  the  test  program  were  operated 
at  a  common  test  condition,  where  P  =  85  psia,  T  =  700  deg  K,  and 
APR  varied  from  50  to  15G,  and  were  found  to  exhibit  significantly 
different  nitric  oxide  emission  characteristics.  At  the  lowest 
AFR  value  their  emission  levels  were  comparable,  but  as  the  AFR 
increased,  the  emissions  from  Combustor  A  decreased  much  more 
rapidly.  This  was  considered  a  convenient  point  to  start  the 
comparisons  between  the  model  predictions  and  the  measured  data, 
as  it  would  clearly  provide"  a  good  test  for  the  model. 

The  first  calculations  were  made,  therefore,  for  this 
test  condition  and  the  range  of  the  mixedness  parameter  Sq  was 
varied  between  0-3  and  0.7  for  both  cases.  The  results, 
compared  in  Figure  16,  show  a  good  agreement  with  experimental 
data  particularly  for  Combustor  B  when  Sq  =  0.7  where  the  predic¬ 
tions  are  approximately  only  20  per  cent  high  over  the  whole  band 
of  AFR  tested.  For  Combustor  A,  the  agreement  is  not  so  successful, 
particularly  at  the  one  point  where  the  air-fuel  ratio  is  approxi¬ 
mately  1A0.  An  acceptable  level  of  agreement  is  obtained  for  the 
case  of  Sq  =  0.7,  but  the  predicted  gradient  of  nitric  oxide  emis¬ 
sions  with  AFR  is  much  lower  than  measured,  whereas  at  Sq  =  0-3, 
the  gradient  is  matched  almost  perfectly  over  the  range  5C  <  AFR^SO, 
but  predicted  values  are  lower  than  the  measured  values.  Before 
discussing  reasons  for  these  discrepancies,  it  is  worthwhile  to 
compare  the  calculated  results  for  these  two  cases  in  more  detail. 

In  Sections  6.1.3  and  6.1.2  it  was  shown  that  the  two  combustors 
under  stuoy  have  significantly  different  design  characteristics. 
Combustor  A  has  a  small  primary  zone  with  a  high  fuel-load  factor, 
hence  low  efficiency,  ar^  operates  very  lean  compared  tc  Combustor  &, 
wnen  operated  at  the  same  AFR  value.  This  iatter  combustor  burns 
relatively  efficiently  in  the  primary  zone,  but  with  a  richer  com¬ 
bustion  mixture.  Calculated  predictions  of  the  internal  details  of 
mass  average  temperature,  equivalence  ratio,  residence  times,  and 
nitric  oxide  levels  are  shown  in  Figure  17  for  AFR  values  o'  about  60. 
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The  comparisons  are  striking.  Combustor  A  is  predicted  Lo  have  a 
top  mass  average  temperature  approximately  275  deg  K  lower  than 
Combustor  B  and  also  a  significantly  shorter  residence  time  to  the 
time  of  freezing  of  the  nitric  oxide  formation  process,  (4  msecs 
versus  8  msecs).  Another  difference  predicted  by  the  model  is  that 
for  Combustor  A  nearly  60  per  cent  of  the  total  NO  emitted  at  the 
exhaust  is  formed  in  the  primary  zone,  whereas  for  Combustor  B, 
only  30  per  rent  has  this  zone  as  its  origin.  This  of  course, 
reflects  the  difference  in  operating  equivalence  ratio  conditions. 

6. 2. 1  Comp  .1  n  with  Combustor  A 

The  disparity  between  the  predicted  and  measured  values  obtained 
for  Combustor  A  and  shown  in  Figure  16  can  be  improved.  It  appears 
that  a  value  of  the  mixedness  parameter  Sq  of  0-3  is  preferable  as  it 
matches  the  slope  of  the  measured  points  very  well.  The  observation 
presented  above,  that  for  this  combustor  the  primary  zone  is  uhe 
active  center  for  NO  production,  suggests  that  a  better  agreement  could 
be  obtained  if  the  primary  zone  fuel-burning  fraction,  /3  ,  were  greater 
than  that  indicated  by  the  combustion  efficiency  correlation.  As 
noted  previously,  the  use  of  this  correlation  tends  to  "under-predi ct" 
(3  .  Therefore,  more  computer  runs  were  conducted,  leaving  all  con¬ 
ditions  constant  except  for  an  increase  in  to  a  value  of  1-15  K/cp  - 
Results  calculated  on  this  basis  show  a  marked  improvement  (see  Fig  18), 
and  now  become  quite  acceptable  for  AFR  values  less  than  about  90.  At 
values  greater  than  this  number,  the  mean  effective  equivalence  ratio 
in  the  primary  zone  soon  drops  lower  than  0-5  (see  Fig  15)  and  due  to 
the  limit  imposed  in  the  model  that  d>  =  2  &>  (see  Volume  3, 

Section  2),  a  region  is  soon  reached  where  no  nitric  oxide  wiil  be 
predicted  (as  St  is  necessary  for  some  mass  to  have  (j>  value  greater 
than  about  0.85  tor  NO  to  be  formed).  Practically,  combustion  becomes 
difficult  to  sustain  at  these  low  values  in  equivalence  ratio  and  the 
difficulty  is  overcome  by  the  use  of  a  pilot  flame  which  operates  from 
its  own  fuel  supply  system.  It  is  clear,  therefore,  that  in  the  experi¬ 
mental  tests,  most  of  the  nitric  oxide  formed  under  these  conditions  wii 
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have  its  origin  at  the  pilot  flame  as  the  pilot  region  is  normally 
designed  to  operate  around  stoichiometric  at  these  high  AFR  values 
for  they  correspond  to  the  idle  operating  conditions  in  aircraft  gas 
turbine  engines.  The  model  clearly  has  its  limitations  for  such  test 
cases,  but  it  is  a  problem  that  can  readily  be  modified  by  accounting 
for  the  pilot  and  main  flows  separately.  As  the  nitric  oxide  levels 
under  discussion  are  only  about  20  ppmV,  .t  is  not  considered  to  be 
a  serious  limitation  at  this  stage. 

Further  calculations  were  made  to  compare  predictions  with  the 
experimental  results  obtained  at  higher  pressure,  P  =  120  psia,  and 
again  the  mixedness  parameter  was  set  equal  to  0.3  and  the  fuel¬ 
burning  fraction  set  equal  to  1.15  l^p  •  The  results  of  these  com¬ 
parisons,  shown  in  Figure  19»  indicate  fair  agreement  between  mea¬ 
sured  and  predicted  emission  rates  at  AFR  values  less  than  90-  However, 
the  model  appears  to  overemphasize  the  effect  of  inlet  air  temperature 
io  that  agreement  is  not  achieved  simultaneously  at  all  temperature 
levels.  it  is  expected  that  more  favorable  agreement  could  be  achieved 
through  further  comparisons  with  variations  in  other  model  parameters. 

6.2.2  Comparison  with  Combustor  B 

a.  Exhaust  Emission  Data.  Figure  16  shows  a  good  agreement  be¬ 
tween  predicted  and  measured  nitric  oxide  emissions  for  Com¬ 
bustor  B  when  the  primary  zone  mixedness  is  assumed  to  equal 
0-7-  As  for  the  case  of  Combustor  A,  it  too  can  be  improved, 
but  for  this  combustor,  as  most  nitric  oxide  is  formed  in  the 
intermediate  zone,  an  improved  agreement  can  best  be  obtained 
by  modifying  the  air  distribution  characteristics-  The  limit 
to  this  modification  clearly  corresponds  to  the  assumption 
that  the  air  entering  the  casing  holes  .nixes  Instantaneously 
with  the  upstream  products.  This  is  not  realistic,  but  does 
serve  to  show  the  sensitivity  of  the  predic  lions  to  the  method 
of  calculation.  Therefore,  ca. eolations  were  made  for  the  cal¬ 
culated  distribution  and  the  1 ii  this  distribu1  on  can  attain 
and  the  results  are  shown  in  Figure  20. 
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Two  observations  can  be  made  about  these  results: 
i.  A  much  better  agreement  is  obtained  at  the  high  AFR 
test  points  than  was  the  case  for  the  previous  com¬ 
bustor  Figure  15  shows  that  for  Combustor  B  at  AFR 
of  120,  close  to  the  highest  value  tested  experimentally, 
primary  zone  effective  equivalence  ratio  is  as  high  as 
0.6  and  thus  it  is  not  surprising  that  the  model  limi-  • 
tations  described  .  6.2.1  do  not  aoply  to  these  test 
conditions. 

ii.  Good  agreement  between  predicted  and  measured  results  is 
obtained  at  one  inlet  temperature  condition  (700  deg  K) , 
but  agreement  is  less  satisfactory  at  other  temperatures.  ■ 
A  different  variation  in  model  parameters  is  required  to 
reduce  the  predicted  sensitivity  to  inlet  temperature, 
b.  Internal  Probe  Test  Results.  The  objective  of  the  internal 
probe  tests  was  to  attempt  to  determine  the  mass  average 
concentrations  of  nitric  oxide  and  carbon  monoxide  as  they 
varied  with  axial  position.  The  basic  data  are  given  in 
Tables  5b  and  5c  and  are  summarized  in  the  table  below,  where 
total  pollutant  mass  flows  are  given  to  be  equal  to  the  product 
of  concentration  and  local  product  flow  rates: 


Measured 

Data 

Predi 

icted  Data 

CO 

CO 

NO 

NO 

NO 

NO 

X 

cone. 

total 

cone- 

total 

cone. 

total 

in 

ppmV 

ppmV  x  M 

3X 

ppmV 

ppmV  x  M 

ax 

ppmV 

ppmV  x  M 

ax 

0.0 

— 

— 

— 

-- 

195 

268 

3-0 

3000+ 

4800+ 

210 

335 

182 

290 

7-0 

425 

1280 

120 

360 

96 

290 

Exit 

120 

710 

58 

345 

49 

290 

Concentration  values  change  with  axial  position  due  to  both  chemi¬ 
cal  reaction  and  dilution  by  the  additional  air  which  enters  the  com¬ 
bustor  liner.  If  there  is  no  chemical  reaction,  then  the  total  mass 
flow  rate  of  the  species  concerned  should  remain  constant  along  the 

.c 

Predicted  data  indicate  that  the  nitric  o\ide  reaction  rate  is 
essentially  zero  at  X  =  2-9  inches. 
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combustor  length.  The  predicted  nitric  oxide  emission  data 
!  indicate:  that  this  occurs  within  the  first  3  inches  of  travel 
■‘from  the' primary  zone.  The  measured  data  for  nitric  oxide  levels 
also,  indicate  the  same  behavior  if  one  makes  allowances  for  the 
1  errors  involved  in.  the  determination  of  the  mass  average  - 
concentration.  This  is  .clearly  not  the  case  for  the  carbon 
monoxide  oxidation  reaction  as  the  experimental  results 
shov  that  there  is  at  least  a  sevenfold  reduction  in  total 
carbon  monoxide  levels  as  the  combustion  products  travel 
from  the  3~inch  station  to  the  exit. 

;  .  =  i 

6. 2. 3  Comparison  at  Low  Power  (idling)  Conditions 

Experimental  tests  were  also  carried  out  on' both  combustors 
at  operating  conditions  which  correspond  to  those  encountered 
in  the  engine  idle  mode.  It  is  known  that  significant  amounts  of 
nitrogen  dioxide  (N02)  can  be  present  in  the  sampled  exhaust  gases 
unden  these  conditions  so  tests  were  made  using  both  wet-chemistry 
techniques: (phenol-disulphonic-acid  method  and  Saltzman)  and  the 
chemiluminescent  analyzer  (see  Section  5-3*3)  for  the  quantity  of 
of  K02  present.  The  results  are  given  in  Figure  21  and  are  also 
compared  with  predictions  made  by  the  model. 

The  following  observations  can  be  made  from  this  figure: 
a.  Total  NOx  levels  as  measured  by  the  Saltzman  technique 
are  significantly  higher  than  the  nitric  oxide  levels 
measured: by  either  the  NDIR  or  chemiluminescent  techniques. 

-b.  Total  NOx  as  measured  by  Saltzman  and  the  chemiluminescent 
analyzer  which  incorporated  an  N02  to  NO  converter  do  not 
agree.  !  The  converter  efficiency  of  the  analyzer  was  checked 
by  passing  known  concentrations  of  N02  through  the  converter 
and  found  to  vary  between  30  and  60  per  cent.  The  discrep¬ 
ancy,  therefore,  appears  to  be  due  to  the  difficulty  in 
conversion  of  the  N02  present, 
c.  Ni r r ; c  oxide  levels  measured  by  the  ND3R  technique  are 
higher  than  those  measured,  wi th  the  chemiluminescent 
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analyzar,  particularly  at  the  high  values  of  AFR.  The 
measurements  made  by  the  two,  however,  tend  to  come  closer 
as  the  AFR  reduces,  that  is,  as  the  combustor  temperatures 
increase.  No  noticeable  discrepancy  occurs  when  the  com¬ 
bustors  operate  at  an  inlet  temperature  higher  than  the 
400  deg  K  that  applies  here.  The  discrepancy  is  'bought 
to  be  due  to  the  high  level  of  unburned  hydrocarbons 
(greater  than  100  ppmV  as  methane)  that  exist  in  the  com¬ 
bustor  exhaust  at  this  condition  as  it  is  known  that  the 
NDIR  technique  does  respond  to  hydrocarbons. 

d.  The  model  predictions  for  nitric  oxide  emissions  agree 
satisfactorily  with  the  nitric  oxide  measurements  made  by 
the  chemi lumi nescent  method. 

6.3  CONCLUSIONS 

The  following  conclusions  are  made  based  upon  the  results  pre¬ 
sented  in  this  section: 

a.  The  model  presented  in  Section  3  adequately  predicts  nitric 
oxide  emission  characteristics  at  high-power  regions  for 
the  two  combustors  tested,  if  one  ignores  operating  condi¬ 
tions  which  produce  nitric  oxide  emissions  of  less  than 

30  ppmV. 

b.  Improved  accuracy  in  prediction  of  nitric  oxide  emissions 
will  depend  as  much  upon  the  accuracy  of  determination  of 
the  primary  zone  mixedness  and  fuel-burning  fraction,  the 
air  distribution  characteristics  within  the  combustor  liner, 
and  the  kinetic  data  as  upon  improvements  to  the  proposed 
flow  models. 

c-  The  model  requires  refinement  if  it  is  to  satisfactorily 

predict  nitric  oxide  emissions  at  low  power  level  conditions 
(basically  defined  by  T  being  less  than  500  deg  K  and  primary 
zone  equivalence  ratio  dropping  below  0.5)*  Consideration 
must  also  be  given  to  nitrogen  dioxide  (NO^)  emissions  which 
constitute  a  substantial  fraction  of  total  NOx  emissions  at 
low  power  levels  (Ref  13)* 
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7-  DESIGN  CRITERIA 

7.1  THE  COMBUSTOR  DESIGN  PROBLEM 

7.1.1*  Performance  Requirements 

In  the  design  of  a  gas  turbine  combustion  chamber,  several  dis¬ 
tinct  performance  requirements  must  be  met.  These  requirements  can 
be  classified  as  follows: 

Aerodynamic  Performance 
a*  Total  pressure  loss. 

b.  Outlet  temperature  profile. 

c.  Wall  temperature  distribution. 

Combustion  Performance 

a.  Combustion  efficiency. 

b.  Combustion  stability. 

c.  Ignitability. 

d.  CO  and  HC  emissions. 

e.  Particulate  emissions  and  smoke. 

f.  NOx  emissions. 

Mechanical  Performance 

a.  Weight. 

b.  Durabil i ty. 
c ■  Cos  t . 

This  method  of  classification  groups  requirements  which  are  most 
closely  related.  Aerodynamic  performance  requirements  are  primarily 
related  to  the  airflow  and  velocity  distributions  in  the  combustor. 
Combustion  performance  requirements  are  related  to  the  detailed  de¬ 
signs  of  the  primary  and  intermediate  zones  and  the  t:uel  injection 
systems.  Mechanical  performance  requirements  are  related  to  the 
mechanical  design  techniques  and  materials  used  in  the  combustor. 
These  performance  groups  are,  of  course,  highly  interrelated. 

In  meeting  the  HOx  emission  control  requ i remen t ,  we  must  be 
concerned  with  its  interactions  with  ^11  the  other  combustion 
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requirements  and  we  find  that  certain  aerodyanmic  requirements  are 
indirectly  affected  by  NOx  emission  control  techniques.  Therefore, 
we  are  concerned  with  those  parameters  which  represent  both  the  com¬ 
bustion  and  aerodynamic  performance  of  the  combustion  chamber  when 
we  design  for  controlled  NOx  emissions.  Mechanical  performance  re¬ 
quirements,  however,  generally  are  not  significantly  affected  by 
combustor  design  modifications  considered  for  NOx  emission  control 
and,  therefore,  will  not  be  included  in  this  discussion. 

7.1.2  Design  Procedures 

The  complete  combustor  design  procedure  is  lengthy  and  complex. 
An  outline  of  the  procedure  is  presented  here  to  identify  the  prin¬ 
cipal  design  parameters  which  the  designer  has  under  his  control. 
This  discussion  applies  to  combustors  of  either  can  or  annular 
configuration.  For  convenience,  terminology  applicable  to  can  con¬ 
figurations  will  be  used. 

tn  designing  a  combustor,  the  following  categories  of  independ¬ 
ent  design  variables  are  involved: 

a.  Combustor  diameter. 

b.  Combustor  length. 

c.  Fuel  nozzle  characteristics. 

d.  Effective  air  hole-  area. 

e.  Conf iguiation  and  location  of  air  holes. 

Th_  ultimate  objective  of  the  design  process  is  the  specif ication  of 
these  quantities.  When  these  design  variables  have  been  specified, 
the  "aerodynamic"  design  of  the  combustor  is  essentially  complete, 
and  the  mechanical  design  process  can  be  initiated- 

To  develop  a  set  of  design  variables  which  w- 1 1  result  in 
sai ;sf aetc.  >  combustor  performance,  the  designer  actually  works  with 
sat  of  intermediate  design  parameters  which  are  more  closely  re¬ 
lated  ^hc  nhysical  processes  occurring  in  the  combustor.  Differ¬ 
ent  combustor  designers  work  with  different  sets  of  intermediate 
parameters,  ^iet  the  set  listed  below  is  representative  of  current 
practice: 
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a-  Pressure  loss  factor. 

b.  Combustion  loading  parameter. 

c.  Combustion  zone  fuel-air  ratio. 

d.  Fuel  nozzle  performance  factors. 

e.  Combustion  zone  length  (or  volume). 

f.  Dilution  zone  length  (or  volume). 

The  pressure  loss  factor  is  essentially  a  combustor  drag  coef¬ 
ficient  and  is  defined  as  follows: 

2  1 

AP  _  AP  f*  "I 

V4  P  L  2RTiJ 


where &P  =  pressure  drop  across  the  combustor 


M 

^ref  P  A«f- 

»  IK 

=  maximum  combustor  casing  cross-sectional  area 


=  pressure  loss  fraction 

n* 

The  pressure  loss  factor  usually  is  assigned  a  value  which  has  been 
shown  from  experience  to  be  consistent  with  satisfactory  combustor 
aerodynamic  performance.  When  the  value  of  this  factor  is  specified, 
the  combustor  diameter  and  pressure  loss  fraction  are  uniquely  related 
for  each  operating  condition. 


The  combustion  loading  parameter  indicates  directly  or  indirectly 
the  heat  release  rate  occurring  in  the  combustion  zone.  Various  forms 
of  this  parameter  have  been  defined  (Ref  15)*  The  form  used  by  NREC 
was  developed  at  Rolls-Royce  Ltd.  (Ref  16)  and  is  defined  as  follows: 


<*}>  (t„A) 
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=  maximum  combustor  casing  diameter 
=  constant  (function  of  combustor  fuel -air  ratio). 
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The  combustion  loading  parameter  is  assigned  a  value  consistent 
with  satisfactory  combustion  performance,  including  combustion  sta¬ 
bility  and  efficiency,  and  ignitability  at  maximum  altitude.  It 
should  be  noted  that  both  the  pressure  loss  factor  and  combustion 
loading  parameter  6 re  related  to  combustor  diameter  so  that  the  com¬ 
bustor  must  be  sized  to  meet  both  criteria- 

Two  fuel -air  ratios  characterize  the  combustion  zone —  that  in 
the  primary  zone  and  that  in  the  intermediate  zone-  Both  quantities 
must  be  specified  in  order  to  fix  the  combustor  airflow  distribution. 

The  performance  of  conventional  fuel  nozzles  can  be  specified 
directly  by  means  of  spray  characteristics —  spray  angle,  velocity, 
and  droplet  size  distribution.  However,  with  the  advent  of  car¬ 
bureting  fuel  injection  systems,  these  characteristics  no  longer 
suffice  since  the  fuel  introduction  process  is  affected  by  the  de¬ 
sign  of  the  combustor  as  well  as  the  nozzle-  A  requirement  exists 
for  a  fuel  injection  performance  parameter,  such  as  the  primary 
zone  mixedness  parameter  (SQ)  defined  in  this  program,  which  indi¬ 
cates  the  distribution  of  fuel  in  the  airflow  and  can  be  related 
to  combustor  and  nozzle  design  variables.  At  present,  the  SQ 
parameter  meets  the  first  requirement,  but  not  the  second.  Thus, 

SQ  cannot  be  specified  quantitatively  at  present.  Hcwever,  mini¬ 
mizing  Sq  js  considered  to  be  an  effective  qualitative  criterion 
for  meeting  certain  combustion  performance  requirements —  notably 
smoke  and  NOx  emission  control. 

The  combustion  zone  length  is  the  sum  of  the  primary  zone  and 
intermediate  zone  lengths.  The  former  is  essentially  fixed  by  the 
combustor  diameter  since  the  reci rculatory  flow  pattern  of  the 
primary  zone  requires  a  "square"  profile.  That  is,  the  primary 
zone  length  must  be  of  the 'same  order  as  the  combustor  radius. 
However,  the  designer  is  free  to  specify  the  length  of  the  inter¬ 
mediate  zone  consistent  with  the  attainment  of  satisfactory  com¬ 
bustion  efficiency,  and  CO  and  HC  emission  rates. 

The  dilution  zone  length  mainly  affects  the  outlet  temperature 
profile.  Thus,  this  factor  is  primarily  related  to  combustor  aero¬ 
dynamic  performance  rather  than  combustion  performance. 


50 


By  specifying  the  values  of  these  intermediate  parameters,  the 
designer  establishes  the  gross  characteristics  of  the  combustor  and 
most  of  the  factors  which  determine  aerodynamic  and  combustion  per¬ 
formance.  The  remaining  tasks  consist  of  detailed  design  of  flow 
passages,  and  specification  of  location  and  configuration  of  cooling 
air  passages. 

The  performance  ar.d  design  parameters  have  been  defined  here  as 
they  are  encountered  in  the  development  of  a  new  combustor.  In  this 
development  problem  which  we  refer  to  as  the  new  combustor  develop¬ 
ment  procedure,  the  performance  factors  are  specified  a  priori  as 
well  as  the  combustor  operating  conditions.  The  design  parameters 
then  are  specified  so  that  the  performance  requirements  will  be 
achieved. 

A  somewhat  different  procedure  is  required  when  an  existing 
combustor  is  to  be  modified,  in  this  development  problem,  which  we 
refer  to  as  the  modified  combustor  development  procedure,  certain 
design  variables,  such  as  length  and  diameter,  are  apt  to  be  fixed 
as  well  as  combustor  operating  conditions.  Under  these  conditions, 
the  constraints  or.  certain  performance  parameters  must  be  relaxed, 
and  an  improvement  in  one  performance  factor  generally  is  accom¬ 
panied  by  a  degradation  of  some  other  aspect  of  combustor  perfor¬ 
mance. 


7-1-3  NOx  Emission  Control 

The  theoretical  analysis  of  NOx  emission  conducted  in  this  pro¬ 
gram  has  indicated  the  dependence  of  NOx  emission  rate  on  engine  pres¬ 
sure  ratio,  combustion  zone  residence  time,  and  fuel  distribution. 

(The  term  "fuel  distribution"  is  used  to  include  both  the  mean  fuel- 
air  ratio  and  the  fuel-air  ratio  distribution  in  the  primary  combus¬ 
tion  zone.)  Of  these  three  factors,  modification  of  fuel  distribu¬ 
tion  is  the  only  effective  approach  to  controlling  NOx  emissions  by 
combustor  design  procedures.  The  engine  pressure  ratio  strongly 
affects  NOx  emission  rate,  but  is  independent  of  combustor  design 
characteristics.  On  the  other  hand,  the  combustion  zone  residence 


time  is  determined  entirely  fay  combustor  design  characteristics, 
but  its  effect  on  NOx  emission  rate  is  relatively  small.  Thus,  it 
is  through  fuel  distribution  that  the  designer  can  most  effectively 
control  NOx  emissions. 

The  fuel  distribution  in  the  combustion  zone  can  be  character¬ 
ized  adequately  by  design  parameters  which  are,  or  could  be,  in¬ 
volved  in  the  combustor  design  process.  These  parameters  are  the 
fuel-air  ratios  in  the  primary  and  intermediate  zones,  and  the  pri¬ 
mary  zone  mixedness  parameter,  Sq.  These  quantities,  therefore, 
are  the  design  parameters  which  must  be  specified  properly  in  order 
to  control  NOx  emissions.  Since  these  parameters  also  affect  other 
combustor  performance  requirements,  they  provide  a  basis  for  relating 
NOx  emission  rate  to  other  combustor  performance  characteristics. 

This  mutual  dependence  upon  combustion  zone  fuel-air  ratio  and  mixed¬ 
ness  underlies  the  development  of  NOx  emission  control  criteria. 

7-2  CONVENTIONAL  COMBUSTORS 

7*2.1  Performance  Characteristics 

It  is  possible  to  depict  the  combustion  performance  characteris¬ 
tics  cf  a  gas  turbine  combustor  in  a  semi -quant i tat ive  manner  by  em¬ 
ploying  certain  design  parameters  defined  in  the  previous  section. 

The  parameters  affecting  combustion  performance  most  strongly  are 
the  combustion  loading  parameter  (© )  and  the  primary  zone  fuel-air, 
or  equivalence,  ratio  (  cpp  ) .  Using  these  parameters,  combustion 
efficiency  can  be  plotted  as  shown  in  Figure  22.  The  figure  con¬ 
stitutes  a  performance  map  for  a  particular  type  of  combustor  and 
indicates  combustici  efficiency  corresponding  to  each  0:  operating 

point.  This  type  of  performance  map  cannot  be  constructed  precisely 
for  a  given  combustor  design  without  the  benefit  of  experimental 
data.  However,  the  map  is  characteristic  of  well-designed  combus¬ 
tors  and  is  useful  for  discussion  of  design  criteria. 

It  can  be  seen  from  Figure  22  that  high  combustion  efficiency 
is  associated  with  high  values  of  ©  and  values  of  Op  near  unity. 
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With  decreasing  ©  and  changes  in  Op  from  this  optimum  value,  com¬ 
bustion  efficiency  is  observed  to  decrease*  The  decrease  is  gradual 
at  first,  but  becomes  precipitous  as  the  efficiency  falls  toward 
low -values  and  the  blow-out  limit.  The  blow-out  limit  of  the  combus¬ 
tor,  as  a  rather  gross  approximation,  may  also  be  considered  to- 
represent  the  combustor  ignition  limit. 

Superimposed  on  each  operating  point,  the  primary  zone  mixed¬ 
ness  parameter,  Sq,  is  shown  as  a  vertical  line.  The  height  of  this 
line  indicates  the  spread  of  equivalence  ratios  over  which  combus¬ 
tion  actually  occurs.  It  can  be  seen  that  if  this  spread  is  great, 
portions  of  the  fuel  are  burning  in  regions  where  combustion  effi¬ 
ciency  is  low  or  where  smoke  is  formed,  even  though  the  mean  equiva¬ 
lence  ratio  is  near  its  optimum  condition. 

The  operating  line  of  a  conventional  (fixed  geometry)  combustor 

l!rf 

is  shown  in  the  figure.  The  locations  of  the  comcustor  operating 
conditions  are  highly  constrained  by  the  requirement  that  the  com¬ 
bustor  operate  over  a  range  of  primary  zone  equivalence  ratios  of 
approximately  2*5:1*  This  range  is  typical  of  the  over-all  fuel-air 
ratio  of  aircraft  turbine  engine  combustors,  and  with  fixed  geometry, 
the  primary  zone  equivalence  ratio  varies  over  a  similar  range.  The 
100  per  cent  power  condition  must  be  located  near  the  =  1  condi¬ 
tion  since  excessive  smoke  is  produced  at  higher  values-  The  value 
of  (|>p  at  100  per  cent  power  cannot  be  reduced  substantially  be¬ 
cause  the  low  power  condition  would  move  out  of  the  stable  combus¬ 
tion  region.  Some  relief  from  this  constraint  is  obtained  by  in¬ 
creasing  ©  through  increases  in  combustor  volume.  However,  the 
dependence  on  0  is  small  in  the  normal  operating  regions  so  that 
any  gain  in  performance  is  small. 

7-2.2  Designing  for  NOx  Control 

Using  the  results  of  the  theoretical  analysis  of  NOx  emissions, 
it  is  possible  to  construct  combustor  performance  maps  depicting  NOx 
emission  rate-  A  different  map  must  be  constructed,  however,  for 
the  new- combustor  and  modif ied-combustor  design  procedures- 

•  he  variable-geometry  operating  line,  also  shown,  is  discussed  in 
Section  7-3-2. 
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In  the  following  discussion,  an  illustrative  example  will  be  given 
of  an  approach  to  NO  emission  control  using  the  modified  combustor 
design  procedure.  Combustor  B  of  the  experimental  test  program  will 
be  used  as  an  example  with  the  assumption  that  the  over-all  dimen¬ 
sions  of  the  combustor  cannot  be  changed.  Approaches  to  NOx  emis¬ 
sion  reauction  will  be  indicated  along  with  the  corresponding  reduc- 
•  tions  in  other  performance  factors. 

It  has  been  shown  earlier  that  NOx  emissions  can  be  controlled 
most  effectively  by  variation  of  the  fuel  distribution  parameters 
(pp  and  Sq.  Thus,  an  NOx  emission  map  for  design  purposes  might  be 
constructed  using  these  parameters.  However,  since  cannot  be 
evaluated  during  the  design  process,  a  different  basis  is  required. 
For  this  purpose,  it  has  been  assumed  here  that  the  primary  zone 
mixedness  factor  Sq  is  related  to  the  pressure  loss  factor  as  fol¬ 
lows: 


c  *AP  ,-l 

°~  Vf 

The  basis  for  this  assumption  is  the  recognition  that  the  pressure 
loss  factor  represents  the  ratio  of  air  jet  energy  to  the  energy  of 
the  main  flow  in  the  combustor.  Thus,  this  factor  is  proportional 
to  the  energy  available  for  "stirring"  the  primary  zone,  and  an  in¬ 
crease  in  stirring  produces  an  increase  in  mixedness,  or  a  decrease 
in  Sq.  The  assumed  relationship  between  Sq  and  &P/qre^  is  a  major 
approximation,  but  no  better  basis  exists  at  present  for  evaluating 
Sq.  Support  for  this  approximation,  based  upon  turbulence  theory, 
is  found  in  Reference  17  in  which  an  approach  to  the  question  of  im¬ 
perfect  mixing  similar  to  the  one  adopted  here  has  been  used,  it 
has  been  shewn  that  the  decay  in  the  factor  Sq  does,  in  fact,  change 
in  proportion  to  the  energy  input  via  air  jets  located  at  the  wall. 
Since  this  discussion  is  intended  only  to  be  illustrative,  this  ap¬ 
proximation  wi  1  1  be  useful.  However,  for  actual  combustor  design 
applications,  a  more  accurate  relationship  between  SQ  and  combustor 
design  variables  will  be  required. 

Exhaust  NOx  concentrations  are  shown  in  Figure  23  for  Combus¬ 
tor  B  operating  at  the  take-off  condition  with  wide  variations  of 


ij>p  and  So»  Using  the  approximate  Sq:  AP/q ref  relationship,  these 
results  have  been  plotted  against  «|i>  andAP/qref  in  Figure  2 b. 

This  figure  can  be  regarded  as  an  NOx  emission  performance  map 
for  this  combustor  since  it  relates  NOx  emissions  to  well-defined 
combustor  design  parameters.  The  map  is  applicable  only  to  the 
modified- combustor  design  procedure  since  it  has  been  constructed 
with  the  assumption  that  over-ail  combustor  geometry  is  fixed. 

The  design  operating  point  for  the  combustor  is  indicated  in 
the  figure.  Also  indicated  is  the  direction  in  which  the  design 
point  might  be  moved  in  order  to  reduce  NOx  emissions-  Since  the  de¬ 
sign  parameters  <pp  andA.P/qre^  are  completely  controlled  by  the  de¬ 
signer,  he  is  free  to  move  the  design  point  about  the  map  as  he 
chooses.  He  must,  however,  be  aware  of  the  consequences  of  his 
modifications. 

It  is  possible  to  relate  three  important  combustor  performance 
parameters—  pressure  loss  f-action,  -ombustion  efficiency,  and  ig- 
nitabllity —  to  the  same  design  parameters —  (jjp  andAP/q^ —  which 
have  been  used  here  to  map  NOx  emissions.  These  relationships  will 
vary  with  combustor  type  and  design  practice,  and  extensive  amounts 
of  performance  data  would  be  required  to  define  them  precisely.  How¬ 
ever,  for  purposes  of  this  discussion,  these  relationships  can  be  de¬ 
fined  in  a  semi-quantitative  manner  from  available  combustor  perfor¬ 
mance  data* 

The  pressure  loss  fraction  is  defined  (Section  7-1*2)  as  a  func¬ 
tion  of AP/q^.  and  combustor  size.  Since  size  is  assumed  fixed  in 
this  example,  the  variation  of  pressure  loss  fraction  is  readily  ob¬ 
served  to  be  proportional  to  variations  inAP/qre^. 

Variations  in  combust:o.i  efficiency  can  be  taken  from  Figure  22. 
With  combustor  size  and  operating  conditions  fixed,  the  loading  parame¬ 
ter  0  remains  fixed,  and  combustion  efficiency  will  vary  with  (pp 
only. 

To  evaluate  variations  in  ignition  performance,  an  ignition 
criterion  developed  by  Lefebvre  (Ref  18)  can  be  used  which  states  that 
constant  ignition  performance  corresponds  to  a  constant  value  of  the 
following  parameter: 
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where  Ej  =  ignition  energy  and  Ure^  is  used  to  represent  the  air  ve¬ 
locity  In  the  combustion  zone.  When  considering  variations  in  (j)p  , 
we  can  substitute  the  inverse  of  this  parameter  for  Ure^  since  is 
changed  by  varying  the  airflow  rate  in  the  combustion  zone.  Then, 
assuming  Ej  is  constant,  the  ignition  parameter  can  be  redefined  as 

X  _  rf.  P°'rT*'r$P  1 
f(Ei) 

The  most  difficult  ignition  requirement  for  an  aircraft  engine  is  re¬ 
ignition  at  high  altitude,  and  the  combustor  inlet  pressure  and  tem¬ 
perature  at  ignition  are  uniquely  related  to  altitude.  Thus,  it  is 
possible  by  the  above  relationship  to  determine  the  variation  in  al¬ 
titude  reignition  capability.  Assuming  that  the  engine  can  be  re¬ 
ignited  at  35,000  ft  altitude  at  present,  reductions  in  this  maxi¬ 
mum  ignition  altitude  can  be  calculated  for  changes  in  and 

^  P/qref * 

The  design  performance  parameter  relationships  discussed  above 
are  presented  graphically  in  Figure  25*  This  figure  is  of  the  same 
format  as  the  NOx  emissions  map  for  Combustor  8.  The  c’  inge  in 
operating  condition,  requi red  to  reduce  NOx  emissions  is  indicated, 
and  the  associated  reductions  in  pressure  loss  fraction,  combustion 
efficiency,  and  reignition  altitude  are  shown.  Also  shown  is  an  es¬ 
timate  of  the  minimum  <f)j»  for  combustor  cooling.  This  estimate  is 
based  on  the  assumptions  that  25  per  cent  of  the  total  airflow  is  re¬ 
quired  for  film  cooling  and  30  per  cent  is  required  for  satisfactory 
dilution  zone  performance,  and  that  the  latter  requirement  decreases 
with  increasing  A?  /Qref 

The  performance  map  in  Figure  25  indicates  that  substantial  per¬ 
formance  losses  would  be  associated  with  reducing  NOx  emissions  by 
variations  in  and  AP/qre^.  The  relationships  shown  in  the  figure 
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are  approximate,  but  they  are  considered  to  be  representative  of  the 
performance  changes  associated  with  the  modif ied-combustor  design 
procedure.  With  the  new-combustor  design  procedure,  where  combustor 
size  is  not  fixed,  the  performance  changes  need  not  be  so  severe. 

The  performance  map  of  Figure  25  indicates  only  the  performance 
losses  at  the  100  per  cent  power  condition.  Performance  losses  also 
would  occur  at  reduced  power  conditions,  particularly  at  ir  e  power. 
Because  of  the  fixed  relationship  between  the  value  of  at  100 

per  cent  and  idle  power  indicated  in  Figure  22,  any  redic  ion  in  dbp 
tends  to  cause  a  rapid  reduction  in  combustion  efficiency  at  idie. 
This  loss  is,  of  course,  accompanied  by  increased  CO  and  HC  emis¬ 
sions.  With  the  modif ied-combustor  design  procedure,  thi:  loss  can 
be  alleviated  somewhat  by  increasing  the  idle  power  level,  but  it  is 
very  likely  that  the  idle  condition  performance  will  be  a  limiting 
factor  if  NO  emissions  are  to  be  controlled  by  reduction  of  4y. 

As  indicated  earlier,  this  discussion  is  intended  to  illustrate 
a  method  by  which  an  existing  combustor  could  be  redesigned  for  low 
NOx  emissions.  If,  instead,  a  combustor  for  a  new  engine  were  to  be 
designed,  variations  in  combustor  geometry  and  fuel  nozzle  character¬ 
istics  could  be  employed  to  minimize  NOx  emissions-  It  is  likely 
that  the  primary  zone  mixedness  parameter,  5q,  is  highly  dependent 
upon  fuel  nozzle  characteristics.  Thus,  a  greater  range  of  design 
conditions  would  be  available  and  the  combustor  performance  losses 
associated  with  NOx  emission  control  probably  would  not  be  as  great 
as  for  a  modified  combustor. 

7-3  ALTERNATIVE  APPROACHES  TO  NOx  EMISSION  CONTROL 

7-3-1  Use  of  Water  Injection 

The  dominant  factor  influencing  the  nitric  oxide  emission  rate 
is  the  peak  temperature  in  the  primary  zone  of  the  rombustor.  This 
temperature  may  be  significantly  reduced  by  the  use  of  water  aodi- 
tion  into  the  combustion  zone,  and  the  water  may  be  added  either 
at  the  compressor  exit  cr  directly  into  the  primary  zone.  The  latter 


approach  is  more  efficient  in  terms  of  the  mass  flow  rate  of  water - 
required  as  approximately  only  one  quarter  of  the  airflow  passes 
through  the  primary  zone. 

With  direct  water  addition  to  t^e  primary  zone,  there  is  a  greater 
possibility  of  a  degradation  of  combustion  performance  than  with  up¬ 
stream  injection.  However,  rec.it  experience  with  direct  water  injec¬ 
tion  {Ref  19)  indicates  that  other  combustion  performance  parameters 
can  be  maintained  while  NOx  emissions  are  reduced.  The  ten.perature 
changes  brought  about  by  water  injection  reduce  the  nitric  oxide 
levels  very  signif ican. !y.  Estimates  have  been  made  of  this  reduc¬ 
tion  for  the  combustors  tested  in  this  program  and  the  results  are 
shown  in  Figures  26  and  27-  The  figures  apply  to  the  case  of  water 
injection  into  the  compressor  exit  and  are  based  uoon  the  experimental 
NOx  emission  data  obtained  by  the  methods  described  earlier.  The  esti¬ 
mates  have  been  made  based  upon  the  following  assumptions: 

a.  That  all  the  water  injected  vaporizes  and  mixes  completely 
with  the  compressor  air  before  entry  into  the  combustor. 

An  enthalpy  balance  then  determines  the  reduction  in  com¬ 
bustor  inlet  temperature  that  results  from  this  action. 

b.  That  the  change  in  nitric  oxide  emission  level  that  results 
from  a  reduction  i.i  temperature  rise  in  the  combustor  (due 

to  the  water),  is  related  to  the  mass  average  flame  tempera-  . 
ture  at  the  primary  zone  exit.  An  enthalpy  balance  accounting 
for  the  presence  of  the  water  vapor  can  determine  the  tempera¬ 
ture  reduction  and  the  corresponding  (higher)  value  of  air- 
fuel  ratio  which  has  the  same  flame  temperature.  Thus,  if  the 
nitric  oxide  versus  air-fuel  ratio  is  known  for  t^e  inlet  tem¬ 
perature  defined  by  assumption  a,  the  nitric  oxide  emission 
level  may  be  estimated  for  a  prescribed  mass  flow  of  water. 

Clearly,  these  calculations  only  produce  estimates  of  the  effect 
of  water  injection  but  these  estimates  will  be  sufficiently  accurate  to 
determine  the  effectiveness  of  the  method.  The  results  are  encouraging 
for  both  combustors  tested  (see  Figs  26  and  27),  particularly  at  the  low¬ 
er  air-fuel  ratios  (less  than  70)  where  a.  water-to-f uel  ratio  of  1.0  gener 
ally  reduces  nitric  oxide  emissions  by  50  per  cent.  The  effect  of  water  i 
jection  is  sho-wn  to  diminish,  however,  as  the  air-fuel  ratio  increases. 


/♦3.2  ■  Variable  Geometry  Combustor  Design  i 

Another. potential ly  effective  approach  to  NO  emission  control 
is  based  on  the  use  of  combustors  with  variable  geometry.  The  pur¬ 
pose  of  the  variability  is  to  maintain  independent  control  of  the 
primary  zone  fuel-air  ratio  as  the  over-all  combustor  fuel-air  ratio 
changes  with  power  level.  This  control  of  {j>p  is  accompl  ished  by 
controlling  the  fraction  of  the  total  airflow  entering  the  primary 
zone.  With  independent  control  of  (|)p  ,  the  combustor  can  be  oper- 
atedso  that  NOx  emission  rates  are  low  over  the  entire  power  range. 
The  most  effective  operating  conditions  for  HOx  control  would  be  at 
lean  fuel-air  ratios  (<j>pO.0).  A  combustor  operating  line  cor¬ 
responding  to  this  condition  is  shown  in  Figure  22.  For  this  ap¬ 
proach  to  be  effective,  a  well-mixed  primary  zone  (low  Sq)  is  re¬ 
quired  so  that  very  little  fuel  is  burned  near  stoichiometric  fuel- 
air  ratios.  An  additional  requirement  is  an  increase  in  combustor 
!size  to  maintain  high  combustion  efficiency  at  the  lower  combustion 
temperatures.-  Because  of  this  requirement,  the  approach  is  not  ap¬ 
plicable  to  existing  engines  where  combustor  size  is  fixed.  But 
variable  geometry  could  be  used  in  future  engines  as  a  means  for 
controlling  NOx  emissions  with  little  or  no  loss  in  other  perfor¬ 
mance  characteristics.  From  the  predicted  effects  of  fuel  distri¬ 
bution  on  NOx. emission  rate,  it  is  estimated  that  reductions  of  50 
to  80  per  cent  can  be  achieved  in  NO  emission  rate  at  high 
power  conditions.  The  approach  also  would  improve  combustion  ef¬ 
ficiency  and  CO  and  HC  emissions  at  low  power,  and  would  improve  ig¬ 
nition  performance.  Thus,  variable  geometry  is  an  attractive  control 
method,  particularly  from  the  standpoint  of  emission  control  effective 
ness. 


7.3*3.  Staged  Combustion 

Another  potential  approach  to  NOx  emission  control  is  based  on 
the  concept  of  staged  combustion  or  staged  fuel  injection.  This  com¬ 
bustor  design  concept  may  take  any  of  several  forms,  but  a  typical 
approach  would  involve  the  following  features: 
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a.  A  small  primary  combustion  zone  operating  at  a  lean  fuel-air 

ratio  (  <j)  0.8)  with  well-mixed  fuel  injection.  The  primary 

zone  would  provide  all  power  required  for  idle  operation,  and 
would  act  as  a  pilot  for  other  combustion  zones  downstream. 

b.  One  or  more  secondary  combustion  zones  into  which  fuel  is 
injected  directly,  again  in  a  well-mixed  fashion.  Fuel  and 
air  would  be  distributed  so  that  fuel-air  ratios  would  be 
lean  throughout  the  combustion  zone. 

This  approach  is  attractive  in  that  it  would  not  require  variable 
geometry,  it  probably  would  require  a  relatively  long  combustor  and 
a  complex  fuel  control  system.  Control  of  NOx  emissions  would  result 
from  the  low  combustion  temperatures  associated  with  lean  fuel-air 
ratios.  The  control  effectiveness  would  increase  with  the  number  of 
combustion  zones  utilized.  However,  tie  fuel  system  complexity  also 
would  increase. 

The  concept  of  staged  combustion  has  not  been  thoroughly  evaluated 
as  an  NOx  emission  control  technique.  Considerable  effort  will  be 
required  to  determine  the  merits  of  this  approach. 


7.4  CONCLUSIONS 

The  following  conclusions  have  been  formulated  concerning  criteria 
for  NOx  emission  control: 

a.  NOx  emission  ratesurelate  to  fundamental  combustor  design 
parameters  which  also  influence  other  major  combustion  perfor¬ 
mance  factors.  This  interrelationship  between  NOx  emissions 
and  other  combustion  performance  factors  establishes  a  sound 
basis  for  quantitative  NOx  emission  control  criteria. 

b.  NOx  emission  rates  from  conventional  combustors  can  be  reduced 
by  modification  of  fuel  distribution  parameters.  However,  NOx 
reduction  may  be  accompanied  by  significant  reductions  in 
other  performance  factors  such  as  combustion  efficiency  and 
reignition  altitude.  These  performance  losses  will  probably 
be  more  severe  in  a  modified- combustor  design  procedure  where 
combustor  size  is  fixed  than  in  a  new-combustor  design  procedure. 
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_.  Water  injection  is  effect iv'  in  reducing  NOx  emission  frn.n 
conventions 5  combustors' .  Emission  rate  reductions  of  50 
per  cent  at  high  power  levels  can  be  obtained  with  water  in¬ 
jection  rates  comparable  to  f>ft]  flow  rates,  increased 
control  may  be  obtained  with  increased  water  injection  rates, 
or  by  localizing  injection  rf  water  into  the  combustion  zone, 
d.  Reductions  ir,  NOx  comparable  to  thnse  achievable  by  the  use 
of  water  injection  can  be  obtained  through  the  development  of 
advanced  combustor  c‘  concepts  such  as  variable  geometry 

and  staged  combustion.  These  concepts  are,  however,  more 
applicable  tc  new  engine  designs  rather  than  existing  designs. 
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TABLE  1 

COMBUSTOR  CONDITIONS  INVESTIGATED  IN 
THE  MODEL  SENSITIVITY  ANALYSIS 


A.  Combustor  Geometry  and  Operating  Conditions 


Voiui.. »  of  Primary  Zone  (V^) 

88.5  in3 

■Length  of  Intermediate  Zone  (X^) 

7-0  in 

Length  of  Intermediate  and  Dilution  Zona 

<W 

12.5  in 

Primary  Zone  Mean  Effective  Equivalence  Ratio 

0.8 

Combustor  inlet  Pressure  (Pj^) 

6.5  atm 

Combustor  Inlet  Temperature  (Tit) 

700  deg  K 

Total  Airflow  to  the  Combustor  (M  ) 

^TOT 

7-56  Ib/sec 

Model  Parameters 

Values 

Tested 

Nominal** 

S  0.1  0.2 

0.3 

0.4 

o 

£*  0.6  0.8 

1.0 

— 

(dMa/dX)o  —  0.2 

0.4 

0.7 

CN*  2  5 

10 

20 

A,*  0.2  0.5 

1.0 

— 

A2*  0.2  0.5 

1.0 

-- 

A3*  0-5  >’.0 

J.5 

— 

Note:  Only  one  variable  was  changed  at  a  time. 


For  definition  of  variables,  see  Section  3*3- 

Nominal  values  for  Aj t  A£,  and  A^  are  the  anticipated  values  based 
upon  the  considerations  given  in  Reference  4.  Nominal  values  for  S 
Cj^  are  estimated.  Nominal  value  for  determined  from  Figure  5»  ° 

and  nominal  value  for  (dM  /dX)  determined  as  described  in  Volume  3, 
Appendix  XVI.  3  c 
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TABLE  2- 

KINETIC  DATA  FOR  THE  NITRIC  OXIDE  REACTION  SCHEME 


-A-here 

k  =  A.t"'  exp(-E./RT) 

i  ) 

T  =  Temperature  in  deg  K 
^  =  1.587  x  10  3  kcal/gm-mole  deg  K 


TABl/i  3 


CHANGES  IN  NITRIC  OXIDE  EM i SSI  ON  LEVELS  PREDICTED 
USING  FOUR  ALTERNATIVE  SETS  OF  KINETIC  DATA 


Kinetic  Data 

1  CODE* 

Reactions** 

1-6 

Reactions** 

1-3 

1 

1.00 

0.98 

2 

2.22 

2.21 

3 

2.05 

2.04 

4 

0.86 

0. 85 

See  Table  2. 

Values  are  quoted  relative  to  1C0DE  =  1  results,  that  is  the  reaction 
data  used  in  Reference  10. 
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TABLE  4a  -  SUMMARY  OF  TEST  DATA 
COMBUSTOR  A 


Combustor  Conditions 


Emission  Levels 


p 

T 

Ma 

AFR 

psia 

deg  K 

deg  R 

Ib/sec 

85 

576 

1037 

5-00 

72 

81 

707 

1273 

3-60 

52 

87 

703 

1265 

4.15 

57 

90 

712 

1282 

3-67 

58 

79 

728 

1310 

3.6I 

64 

91 

705 

1269 

4.13 

70 

82 

687 

1237 

3-43 

71 

84 

700 

1260 

3.51 

76 

88 

698 

1256 

4.00 

79 

89 

713 

1283 

4.23 

83 

75 

683 

1229 

3-85 

143 

89 

8.6 

1469 

3.19 

63 

83 

688 

1238 

2.74 

59 

85 

711 

1280 

5.03 

61 

119 

579 

1042 

6.2? 

56 

1.48 

598 

1076 

6.31 

67 

118 

581 

1046 

5.61 

71 

123 

578 

1040 

5-89 

75 

116 

590 

1062 

6.63 

82 

114 

595 

1071 

5-67 

95 

112 

603 

1085 

6.16 

153 

123 

675 

1215 

5-72 

54 

119 

688 

1238 

5-57 

59 

113 

677 

1218 

5-02 

64 

130 

692 

1245 

5-64 

67 

131 

698 

1256 

5-75 

72 

110 

691 

1243 

4.97 

74 

108 

695 

1251 

5-27 

10^ 

102 

695 

1251 

5-06 

116 

MaT/P 

NO 

ppmV 

CO 

ppmV 

HC 

ppmV 

(M) 

34 

46 

480 

— 

1 

31 

184 

— 

— 

1 

34 

137 

245 

— 

1 

29 

160 

206 

0.8 

1 

33 

134 

250 

— 

1 

32 

102 

275 

— 

1 

29 

101 

340 

— 

1 

29 

80 

— 

— 

1 

32 

58 

— 

— 

1 

34 

34 

— 

— 

1 

35 

21 

245 

1 

29 

192 

170 

— 

1 

23 

136 

125 

— 

1 

42 

110 

190 

0.6 

30 

111 

255 

— 

1 

32 

83 

— 

— 

1 

28 

58 

290 

1.3 

1 

28 

56 

233 

1.3 

1 

34 

56 

380 

— 

1 

30 

32 

300 

3. 1 

1 

33 

19 

300 

3.1 

1 

31 

264 

205 

— 

1 

32 

228 

195 

— 

1 

30 

122 

182 

— 

1 

30 

122 

255 

2.0 

1 

31 

104 

265 

2.0 

1 

31 

88 

210 

— 

1 

24 

53 

— 

0.9 

1 

34 

44 

137 

0.9 

1 

Fuel:  Aviation  Kerosene 


Calibrated  as  hexane 

Method  of  NO  measurement:  1  =  ND1R 

2  =  Chemiluminescent 
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TABLE  4b  -  SUMMARY  OF  TEST  DATA 
COMBUSTOR  A 


Combustor  Conditions  Emission  Levels 


p 

T 

Ma 

AFR 

MaT/P 

NO 

Cu 

h'C 

(M) 

psia 

deg  K 

deg  R 

lb/ sec 

ppmV 

ppmV 

ppmV 

123 

820 

1476 

4.09 

53 

27 

376 

120 

1.2 

1 

123 

852 

1533 

4.04 

53 

28 

340 

— 

1.3 

1 

124 

845 

1521 

4.25 

66 

29 

269 

167 

1.2 

1 

119 

874 

1573 

4.15 

66 

30 

254 

174 

1.3 

1 

*  *  — 

•  •  *3 

855 

1539 

4.11 

68 

30 

213 

— 

1.3 

1 

117 

852 

1533 

4.25 

89 

31 

155 

133 

0.8 

1 

36 

402 

724 

0-93 

78 

10 

33 

760 

Sat. 

1 

42 

401 

722 

0.91 

44 

9 

73 

1420 

Sat. 

1 

Fuel:  Aviation  Kerosene 


Calibrated  as  hexane 
Method  of  NO  measurement: 


1  =  NDIR 

2  =  Chemiluminescent 
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TABLE  5a  -  SUMMARY  OF  TEST  DATA 
COMBUSTOR  b 


Combustor  Conditions  Emission  Levels 


p 

T 

Ma 

AFR 

MaT/P 

NO 

CO 

HC 

(M) 

psia 

deg  K 

Deg  R 

lb/sec 

ppmV 

ppmV 

ppmV 

•WU 

89-7 

5<?3 

1049 

6. 10 

70 

40 

81 

120 

1.0 

1 

89.7 

583 

1049 

6. 10 

75 

40 

72 

— 

— 

2 

87-2 

576 

1036 

6. 10 

86 

40 

69 

67 

1.0 

1 

87-2 

576 

1036 

6. 10 

87 

40 

65 

— 

— 

2 

86.7 

588 

1058 

6.16 

92 

42 

69 

113 

2.0 

1 

86.7 

588 

1058 

6.16 

92 

42 

63 

— 

— 

2 

87-6 

593 

1067 

5-80 

113 

39 

57 

93 

— 

1 

88 

588 

1058 

5-97 

126 

40 

34 

134 

1,2 

91 

683 

1229 

5- 16 

56 

39 

153 

87 

— 

1 

91 

683 

1229 

5- 16 

56 

39 

148 

— 

2 

86 

687 

1236 

4.99 

62 

40 

152 

100 

— 

1 

86 

687 

1236 

4-99 

62 

40 

149 

— 

— 

2 

83 

703 

1265 

5.04 

75 

43 

148 

87 

— 

1 

83 

703 

1265 

5-04 

75 

43 

144 

— 

2 

82 

703 

1265 

5-09 

96 

44 

115 

80 

— 

1 

82 

703 

1265 

5  J9 

96 

44 

113 

— 

2 

80 

708 

1274 

5.31 

129 

47 

76 

SC 

— 

1 

80 

708 

1274 

5.31 

129 

47 

76 

—— 

2 

82 

797 

1434 

4.61 

66 

45 

235 

100 

1.0 

1 

82 

797 

1434 

4.61 

71 

45 

215 

— 

— 

2 

85 

797 

1434 

5-03 

93 

47 

I69 

67 

1-5 

1 

85 

797 

1434 

5-03 

100 

47 

150 

— 

— 

2 

82 

827 

1488 

4.77 

107 

48 

169 

67 

3-5 

1 

82 

827 

1488 

4.77 

115 

48 

144 

— 

— 

2 

37 

393 

707 

3-66 

127 

30 

14 

652 

— 

1 

37 

393 

707 

3-66 

127 

39 

3 

— 

— 

2 

28 

392 

705 

3-62 

76 

51 

20 

500 

— 

1 

28 

392 

705 

3.62 

76 

51 

16 

— 

-- 

2 

Fuel:  Aviation  Kerosene 


Calibrated  as  hexane 

Method  of  NO  measurement:  1  =  ND1R 

2  =  Chemiluminescent 
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TABLE  5b  -  SUMMARY  OF  TEST  DATA 
COMBUSTOR  B 


Interna]  Probe  Measurement,  X  =  3  inches 
Combustor  Conditions  Emission  Levels 


r 

i  nches 

p 

ps  ia 

deg  K 

T 

deg  R 

Ma 

lb/sec 

AFR 

MaT/p 

NO 

ppmV 

CO 

ppmV 

(M) 

Exit 

71 

595 

1071 

5-75 

85 

48 

68 

98 

2 

Wall 

71 

600 

1080 

5-86 

88 

50 

133 

1750 

2 

3-0 

71 

59-3 

1067 

5-81 

88 

49 

143 

2050 

2 

2.5 

72 

598 

1076 

5.91 

89 

49 

150 

2400 

2 

2.0 

71 

578 

1040 

5.88 

89 

48 

156 

3000+ 

2 

1-5 

70 

59^ 

1069 

5.80 

88 

**9 

216 

+ 

2 

Exit 

63 

588 

1C58 

5-83 

88 

50 

56 

— 

2 

Exit 

66 

600 

1080 

5*82 

88 

53 

60 

150 

2 

2.0 

68 

605 

1089 

5.99 

90 

53 

162 

2200 

+• 

1.5 

68 

583 

1049 

5.9** 

89 

51 

205 

+ 

2 

1.0 

67 

581 

1045 

5.88 

88 

51 

278 

+ 

2 

0.5 

68 

583 

1043 

5.93 

88 

51 

305 

+ 

2 

0.6 

67 

581 

1045 

5-78 

86 

50 

340 

+ 

2 

Er  i  t 

68 

595 

1071 

5.93 

88 

52 

66 

— 

.2 

Approximate  values  of  mass  average  concentrations  of  pollutant  species 
calculated  to  be:  NO  =  211,  CO  =  7000+ 


Fuel:  Aviation  Kerosene 


Method  of  NO  measurement: 
^  Distance  from  center  line 
Assume  equilibrium  values 


1  =  ND1R,  2  =  Chen i luminescent 


SUMMARY  OF  TEST  DATA 
COMBUSTOR  B 


Internal  Probe  Measurements,  X  =  7  inches 


Combustor  Conditions 


Emission  Levels 


r 

P 

T 

M 

AFR 

M  T/P 

NO 

CO 

(M) 

i nches 

psia 

deg  K 

deg  R 

lb/sec 

O 

ppmV 

ppmV 

Exit 

73 

599 

1078 

5-94 

87 

49 

59 

120 

2 

Wall 

71 

593 

1067 

5.86 

86 

49 

33 

340 

2 

3-0 

71 

595 

1071 

5-86 

87 

49 

56 

400 

2 

2.5 

72 

575 

1035 

5-98 

89 

48 

115 

620 

2 

Exit 

70 

596 

1072 

5-85 

87 

50 

54 

146 

2 

2.0 

71 

581 

1045 

5.95 

89 

49 

125 

486 

2 

1.5 

71 

584 

1051 

5-98 

90 

49 

112 

340 

2 

Exi  t 

72 

594 

1069 

5.98 

90 

49 

54 

142 

2 

1.0 

71 

592 

1065 

5-95 

90 

50 

122 

333 

2 

0.5 

72 

593 

1067 

6.03 

91 

50 

132 

340 

2 

0.0 

70 

591 

1063 

6.03 

92 

51 

— 

433 

2 

Approximate  value  of  mass  average  concentrations  of  pollutant  species 
calculated  to  be:  NO  =  120  rpm,  CO  =  425  rpm 

Fuel:  Aviation  Kerosene 


Calibrated  as  hexane 

Method  of  NO  measurement:  l  =  NDIR,  2  =  Chemiluminescent 
Distance  from  center  line 
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PRIMARY  INTERMEDIATE  DILUTION 


(a)  CROSS  SECTION  of  COMBUSTOR  LINER 
SHOWING  MEAN  FLOW  PATTERN 


(bl  DISTRIBUTION  of  MIXTURE  RATIO  F, 
ALONG  LINER  LENGTH 


F  AT  POSITION  A 


(c)  DISTRIBUTION  of  MASS  AS  A  FUNCTION 
of  MIXTURE  RATIC  and  POSITION 


Figure  1  -  ASSUMED  FLOW  AND  DISTRIBUTION  PATTERNS 
THROUGHOUT  COMBUSTOR  LINER 


ANO/NOnnm ,  PERCENT 


(o).  EFFECT  of  /? 


(b).  EFFECT  of  dM„/dX 


Figure  3 -EFFECT  OF  PRIMARY  ZONE  FUEL-BURNING  FRACTION 
AND  DILUTION  AIR  MIXING  RATE 
UPON  PREDICTED  NITRIC  OXIDE  EMISSIONS 


NO  lb/ 1000  lb  FUEL  BURNED 


MODE 

nom(msecs) 

TAKE  OFF 

2.40 

CRUISE 

2.55 

APPROACH 

2.63 

IDLE 

3.  13 

*  ALL  IDLE  VALUES  ARE 
LESS  THAN  0.01  Ib/lOOOlb 
FUEL  BURNED. 


6 

5 


0.6  0.8  1.0  1.2  1.4 

TP  ^TP  nominal 


Figure  5-  EFFECT  OF  PRIMARY  ZONE 
RESIDENCE  TIME  ON  NGX  EMISSION  PREDICTIONS 
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8 


EPR 

RADIUS  tin}  1 

12=1 

r.33 

!6.?=e 

3.25 

25:1 

2.65 

*  ALL  IDLE  VALUES  ARE 
LESS  THAN  0.01  Ib/lOOOlb 
R-€t  BURNED. 
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IDLE* 
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16  20  22  24  26 


14  16 

EPR 

(RATED  PRESSURE  RATIO  AT  FULL  POWER} 


Figura  7  -  FFFECT  OF  ENGINE  PRESSURE 
RATIO  UPON  NO.  EMISSION  PREDICTIONS 
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DATA  FROM  REFERENCE  13 
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EPR 


(Rated  Pressure  Ratio  at  Full  Power ) 

Figure  e  -  MEASURED  EFFECT  OF  ENGINE 
PRESSURE  RATIO  UPON  N0X  EMISSIONS 


FUEL  LOAD  FACTOR 
Mf  /Vp  P2  -lb /sec  ft 3  Atm2 

REPRODUCED  FROM  REFERENCE  15 

Figure  9  -  PRIMARY  ZONE 
COMBUSTION  EFFICIENCY  CORRELATION 


Figure  io- COMBUSTOR  TEST  ARRANGEMENT 


84 


Figure  12  ■  THE  TEST  CONTROL  ROOM 


s& 


BLEED  TO 

NO  CHEMILUMINESCENCE  ANALYSER 


Figure  <3 -SCHEMATIC  DIAGRAM  of  EQUIPMENT  EMPLOYED  IN 
EMISSIONS  SAMPLING  and  ANALYSIS 


Figure  u-A'R  DISTRIBUTION  CHARACTERISTICS 


PRIMARY  ZONE  EQUIVALENCE  RATIO 


8‘ 


Figure  15- PRIMARY  ZONE  OPERATING  CHARACTERISTICS 


NITRIC  OXIDE  EMISSION,  ppmV 


OPERATING  (  P  =  85  psia 
CONDITIONS  \  T  =  700  deg  K 

COMBUSTOR  CXP.  THEORY 

A  — £r I  1 

B  — O”  —  1-1 


40  60  80  100  120  140  I6C 

AIR-FUEL  RATIO,  AF° 

S0  -MIXEDNESS  IN  PRIMARY  ZONE 


Figure  16 -PREDICTED  versus  MEASURED  NITRIC  OXIDE  EMISSIONS 


NITRIC  OXIDE ,  ppmV 
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OPERATING  P  =  83  p*io 
CONDITIONS  T  =  710  deg  K 


A  MEASURED 

- PREDICTED  0-  77cp  (FROM  FIG.  9) 

-  PREDICTED  0  =  1. 15  if 


Figure  *8-  PREDICTED  versus  MEASURED  NITRIC  OXIDE  EMISSIONS 

(COMBUSTOR  A,  P  =  85 ) 


NITRIC  OXIDE ,  ppmV 


OPERATING  I  p 
CONDITIONS  J  V 


120  psia 
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MEASURED  PREDICTED  (So  =  0.3,  £  =  I-I5  i)  ) 

T  =  600  □  P 

700  A  - 

850  O  - 


AIR-FUEL  RATIO ,  AFR 
S0  -NIXEDNESS  IN  PRIMARY  ZONE 


igure  is-PREOICTED  versus  MEASURED  NITRIC  OXIDE  EMISSIONS 

(COMBUSTOR  A,  P=  120) 


.1 


1.0  !0 
COMBUSTION  LOADING  PARAMETER* 10)  ~,if08 


0= 


p  T75  0.75 

•  Aref  D ref 


e,p(- tJ 


*  TYPICAL  VALUES  PLOTTED  IN  lbm,  ibj.ft,  sec,  deg  R  UNITS. 


Figure  22  -  COMBUSTION  PERFORMANCE  MAP  FOR 
FIXED  and  VARIABLE  GEOMETRY  COMBUSTORS 


9' 

P  =  17  Aim 
T  =  700  deg  K 
AF*  =  68 


SIMULATED  TAKE-OFF 


PRIMARY  ZONE  EQUIVALENT  RATIO 


S^Sj£3*Se2; 
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?  =  17  Atm  j 

T  =  700  deg  K  j  SIMULATED  TAKE-OFF 
AFR  =  68  J 


PRESSURE  LOSS  FACTOR  AP/g 

Figure  25  -  PERFORMANCE  MAP  FOR  MODIFIED  -  COMBUSTOR 
DESIGN  PROCEDURE  (COMBUSTOR  B) 


000  lb  FUEL 


COfoBUSTOS  A. 

INLET  T  s  850  deg  K 
ps  !20  psia 


WATER  INJECTED  INTO  COMPRESSOR  AIR 


MASS  WATER 
mass  FUEL 


NITRIC  ..OXIDE  EMISSIONS ,  lb/ 1000  lb  ENEL 


COMBUSTOR  B 
INLET  T  =  800  deg  K 


AIR-FUEL  RATIO 

-  t 

1  .  i 

figure  2? -ESTIMATED  EFFECT  OF  WATER  INJECTION  (COMBUSTOR  B) 
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1. 


NOMENCLATURE 


Symbols 


A 

E 

F 

k 

M 

P 

R 

S 

o 

t 

T 

T 

m 

V 


r 

e 


$ 

p 


4>. 


max 


icp 


•  Subscripts 

a 

f 

i 

o 


Description 

Preexponential  factor  in  Arthenus  equation 
Total  cross-sectional  area  of  combustor  liner 
Activiation  energy 
Mixture  ratio  (fuel  mass/total  mass) 

Reaction  rate  constant 

Mass  flow  rate  in  combustor  liner 

Pressure 

Radius  of  combustor  liner 

Degree  of  mixedness  in  primary  zone  («=  CT^/Fp) 

Time 

Temperature 

Maximum  cycle  temperature 
Volume 

Axial  distance 

Length  of  intermediate  zone 

Distance  from  primary  zone  exit  to  the  liner  exit 
Fraction  of  fuel  burned  in  primary  zone 

WVWItwOl  V  VVV  IOvi 

Dens i ty 

Residence  Time 

Combustion  efficiency  parameter 
Equivalence  ratio  in  primary  zone  **  1) 

Effective  equivalence  ratio  in  primary  zone  (=$^p) 

Maximum  equivalence  ratio 

Primary  zone  combustion  efficiency 

Description 

Ai  r 

Total  fuel 

i‘ th  element  in  series 

Condition  at  entrance  to  intermediate  zone 


10) 


Subscripts 

P 

x 

in 


Description 

Condition  in  primary  zone 

Axial  position 

Condition  at  combustor  inlet 
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